

LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
KIRTLAND AFB, N MEX 


TRANSIENT THERMODYNAMIC ANALYSIS 
OF A FUEL-CELL SYSTEM 

by William Emile Simon 

Manned Spacecraft Center 
Houston, Texas 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • JUNE 1968 



TECH LIBRARY KAFB, NM 



tech library kafb, nm 



NASA TN D-4601 


TRANSIENT THERMODYNAMIC ANALYSIS 
OF A FUEL-CELL SYSTEM 

By William Emile Simon 

Manned Spacecraft Center 
Houston, Texas 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 




ABSTRACT 


This report presents a transient thermodynamic 
analysis of a fuel-cell system exclusive of cryogenic 
supply system and external heat dissipation system. 
The results of a simplified analysis and of a more 
general analysis are compared, and the results of the 
more general analysis are compared with perform- 
ance test data. The analysis is oriented toward a 
particular application, but the methods employed are 
sufficiently general to allow application of these con- 
cepts to most fluid-loop systems involving heat and 
mass transfer. 


CONTENTS 


Section Page 

SUMMARY 1 

INTRODUCTION 1 

SYMBOLS 3 

OPERATING PRINCIPLES AND CHARACTERISTICS OF THE FUEL CELL ... 9 

DESCRIPTION OF FUEL-CELL SYSTEM SELECTED FOR ANALYSIS 17 

GENERAL ANALYTICAL APPROACH 22 

Type of Transient Analysis 23 

Current Variations 25 

Fuel-Cell Environment Temperature 25 

Thermal Mass 25 

Chemical Equilibrium 25 

Stoichiometry 26 

Perfect-Gas Law 26 

Condenser- Exit- Temperature Variations '*. . . 26 

Pressure Conditions 26 

Primary-Loop Volumetric Flow Rate And Efficiency of Water 

Separation 27 

Performance Degradation 27 

Reactant Inlet Temperature 27 

Homogeneity of Hydrogen and Water-Vapor Mixture 27 

DEVELOPMENT OF GOVERNING EQUATIONS 27 

SIMPLIFIED TRANSIENT ANALYSIS 39 



Section 


Page 

GENERAL TRANSIENT ANALYSIS 41 

Case I 49 

Case II 53 

Case III 57 

Case IV 58 

Case V 59 

COMPARISON OF ANALYTICAL RESULTS 60 

CONCLUDING REMARKS 63 

APPENDIX A - DEVELOPMENT OF EMPIRICAL RELATION FOR FUEL-CELL 
POWER PLANT OUTPUT VOLTAGE AS A FUNCTION OF 
CURRENT, TEMPERATURE, AND ELECTROLYTE 
CONCENTRATION 66 

APPENDIX B - THE COMPUTER PROGRAM FOR THE GENERAL 

ANALYSIS 71 

REFERENCES 81 


IV 


TABLES 


Table Page 

I CASES CONSIDERED IN THE DEVELOPMENT OF THE TRANSIENT 
EQUATIONS FOR DIFFERENT CURRENT AND TEMPERATURE 
REGIONS OF OPERATION 49 


v 



FIGURES 


Figure Page 

1 Fuel-cell heat and mass balance 10 

2 Energy diagram for the fuel cell 11 

3 Hydrogen-oxygen fuel-cell performance based on equilibrium 

conditions 12 

4 Fuel-cell chemical balance 14 

5 Fuel- cell electrochemical flow schematic 15 

6 Reaction mechanisms at hydrogen and oxygen electrodes 

(a) Hydrogen electrode 16 

(b) Oxygen electrode 16 

7 Apollo fuel- cell module 11 

8 Energy conversion section 18 

9 Accessory section — Bay I 18 

10 Accessory section — Bay H 19 

11 Accessory section — Bay IH 19 

12 Single-cell diagram illustrating cell construction 20 

13 Schematic diagram of fuel-cell system 21 

14 Control-volume approach to the analysis 23 

15 Transient fuel-cell performance 24 

16 Transient fuel- cell performance (step current increase) 24 

17 Transient fuel- cell performance (step current decrease) 24 

18 Fluid- flow diagram and heat balance for fuel- cell stack 

(a) Flow diagram 29 

(b) Heat balance 29 



Page 


Figure 

19 Specific heats at constant pressure for oxygen, hydrogen, and 

water vapor 32 

20 Primary-bypass- valve characteristics 46 

21 Comparison of results of the simplified analysis and the 

general analysis 60 

22 Comparison of steady-state voltage-current characteristics 61 

23 Comparison of results of general analysis with vacuum- chamber 

test data 

(a) 159.9° F ^ T. <166.3° F 62 

(b) 157.0° F < T, < 167. 1° F 62 

(c) 158.0° F < T^ e < 165. 5° F 62 

(d) 159. 5°F< T c y e < 164. 5° F 62 

A-l Comparison of polarization curves of analytical model with those 

of statistical test data 66 

A-2 Variation of slope and Y-intercept of constant-temperature 

voltage lines with temperature 67 

A- 3 Variation of electrolyte concentration with temperature and 

partial pressure 69 

A-4 Variation of Y-intercept of constant-temperature concentration 

lines with temperature 69 


vii 



TRANSIENT THERMODYNAMIC ANALYSIS 

OF A FUEL-CELL SYSTEM 

By William Emile Simon 
Manned Spacecraft Center 


SUMMARY 


This report is concerned with the development of an analytical model capable of 
accurately predicting the transient-temperature and voltage-response characteristics 
of a fuel-cell system. The Apollo fuel-cell system was selected for analysis. For a 
given space mission, a complete time history of fuel-cell temperature and voltage must 
be obtained before launch to determine (1) whether there are sufficient reactants aboard 
the spacecraft available for fuel-cell consumption to complete the objectives of the mis- 
sion and (2) whether fuel-cell voltage is always within the specification voltage band of 
29 ± 2 volts. The amount of reactants consumed depends on the efficiency of the fuel 
cell, which is a function of operating temperature. 

To obtain the necessary relations for the variation of fuel-cell temperature and 
voltage with time, a control volume was set up around the fuel-cell stack, which is the 
heart of the system. Heat and mass balances were performed, and appropriate bound- 
ary and initial conditions were applied. Because of the complexity of the system, it 
was necessary to introduce a set of simplifying assumptions to hold the number of 
varying parameters to a minimum. Some empirical (test) data were also needed. 

After the transient equations were developed for all regions of normal operation, 
these relationships were used to formulate a computer program which accepts a given 
mission profile and calculates the temperature and voltage variations of the fuel cell 
with time. 

From a comparison of analytical results with test data obtained on a qualified 
Apollo fuel-cell power plant in a thermal vacuum environment, it was concluded that the 
analysis presents an adequate mathematical model for predicting fuel-cell performance, 
subject to the limitations imposed on the study. 


INTRODUCTION 


Analytical techniques directed toward the prediction of system performance are 
important in the design of all systems. The prediction of the operating characteristics 
of the system is based on the physical geometries of the system components and on the 


behavior of the working media. When various forms of energy are converted, noniso- 
thermal systems may result. Thus, the performance analysis may be significantly 
complicated because (1) the system performance is generally temperature dependent, 
and (2) heat transfer mechanisms may be present. Generally, a thermal analysis ex- 
tends beyond the system to include the influence of the surroundings. If the system is 
in thermodynamic equilibrium, it is in the steady- state condition. If any of the vari- 
ables of the system are changing with time, the system is in the transient state. 

Accurate analysis of system performance is particularly critical for design con- 
siderations for systems such as those involved in space vehicles. This report is con- 
cerned with the development of an analysis to describe the performance of a fuel-cell 
system such as the one used in the Apollo spacecraft. The fuel-cell system aboard the 
Apollo spacecraft accepts hydrogen and oxygen from a cryogenic storage and supply 
system and converts these reactants electrochemically into waste heat, water, and 
electrical energy. A thermal- control system removes the waste heat for dissipation 
into space by spacecraft radiators. 

The variables which govern the performance of a multicomponent fuel-cell system 
are numerous. Analytical treatments rely to a certain extent on empirical (test) data 
and are limited by necessary simplifying assumptions. The most important parameter 
governing the performance and life capability of this type of fuel cell is temperature. 
Although performance increases with increasing temperature, fuel-cell life is inversely 
proportional to temperature. The voltage output of the system is highly dependent on 
fuel-cell temperature at any given load. Since the allowable range of voltage variation 
is very sensitive (only ±2 volts about a nominal 29 volts), an accurate thermal analysis 
is necessary to be able to predict the performance of a fuel-cell system for a given 
space mission. 

Thermodynamic analyses of fuel-cell systems are largely limited to steady-state 
treatments and can generally be classified as in-house reports. A transient analysis 
must determine the proper relationships for the variation of fuel-cell temperature and 
voltage with time after a given step load change for all regions of operation. Inputs to 
an analysis of this type must include (1) load profiles (current as a function of time), 

(2) fuel-cell coolant inlet-temperature profiles, (3) reactant inlet-temperature profiles, 
and (4) variations in the environment temperature of the fuel-cell system. The desired 
outputs of the analysis are the transient temperature and the voltage response of the 
fuel cell in terms of such variables as electrolyte concentration. 

This report presents a transient thermodynamic analysis of a fuel-cell system 
exclusive of cryogenic supply system and external heat dissipation system. The results 
of a simplified analysis and of a more general analysis are compared, and the results 
of the more general analysis are compared with performance test data. The analysis is 
oriented toward a particular application, but the methods employed are sufficiently gen- 
eral to allow application of these concepts to most fluid-loop systems involving heat and 
mass transfer. 
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SYMBOLS 


A, a, B, b, .C, c 
BP1, BP2 

b l 

C 1 

CAHI, CBHI, CDHI, CEHE 


CALO, CBLO, CDLO, CELO 


E 


E, 


E f 


E 


E 


cs 


E 


rev 


E 


ss 


E 


e 

G 

H 


collection of terms used in the analysis for simplification 

curve fit constants for primary-bypass-valve character- 
istics 

Y-intercept 

constant used in equations (135), (137), and (138) 

curve fit constants for fuel-cell performance character- 
istics in the high- current region, 5 A ^ I s 55 A 

curve fit constants for fuel-cell performance character- 
istics in the low- current region, 0sI<5A 

specific heat at constant pressure, Btu/lb m -°R 

electromotive force of the fuel cell or module output 
voltage, V 

initial voltage, V 

second-stage or final voltage, V 

actual electromotive force, V 

cold- side effectiveness of primary regenerator 

reversible electromotive force, V 

steady- state module voltage, V 

voltage variation with current and temperature 
(appendix A), V 

base of natural logarithms 

free energy (Gibbs function), Btu 

heat content (enthalpy of formation), Btu/° R 
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HF 

h 



heat of formation of water, Btu/lb -H 0 

’ 'm2 

specific enthalpy, Btu/lb m -°R 

latent heat of water condensed, Btu/lb -H o 0 

’ m2 

module current, A 
initial current, A 

final current after a step load change, A 


2 2 
I R heating load, A 

I. , total current, A 

tot 

Kj electrochemical constant defined by Faraday's law, 

A-hr/lb -H 0 cons 
m 2 


K correction factor for the electrolyte constant of the fuel 

c cell, V/A percent HgO above 27 percent 

M thermal mass of the fuel-cell stack, Btu/°F 

MW molecular weight, lb m /lb-mole 


m 


mass, lb 

m 


m 


mass flow rate, lb /hr 
’ m 


m^ slope 

n motor speed, rpm 

P gross electrical power output of the fuel-cell module, W 

2 

P . I R loss in motor, W 

elec 

Pj w combination of friction and windage losses of power, W 

P combination of hysteresis and eddy-current losses of 

he power, W 


4 



power dissipated in motor because of losses, W 


P L 

P shaft 


P 

P H 

P w 

Q 


QA, QB 

q a 

Q 

gen 

Q htr 

Ql 

Q 

recir 

Q 

^s 

^stack 

q 

R 

R 

S 


shaft output power of motor, W 

steady- state power level, W 

absolute (or total) pressure, psia 
partial pressure of hydrogen, psia 

partial pressure of water vapor, psia 

heat generated because of the inefficiency of the fuel cell, 
adjusted for the difference in temperature between 
incoming reactants and outgoing products, Btu/hr 

curve fit constants for stack heat loss 

heat produced in an actual fuel cell, Btu 

product of hydrogen- consumption flow rate and lower 
heating value, Btu/hr at 77° F and 1 atm 

in-line heater power, Btu/hr 

portion of waste heat generated in actual cell because of 
losses, Btu/hr 

heat removed by the recirculating stream through the 
stack, Btu/hr 

heat stored in stack, Btu/hr 
stack heat loss, Btu/hr 

heat quantity (also, in some parts of the analysis, a 
collection of terms used to form a discriminant), Btu 

gas constant, ft-lb./lb -°R 
f m 

universal gas constant, ft-lb f /lb m -mole-°R 
entropy, Btu/°R 
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w 

X 


specific entropy, Btu/lb m -°R 

fuel-cell operating temperature, °F or °R 
reference temperature, °F or °R 

initial fuel-cell temperature, °F or °R 

bypass temperature, high 

bypass temperature, low 

condenser exit temperature, °F 

upper limit of in-line (low-power) heater temperature 
dead band, °F 

lower limit of in-line (low-power) heater temperature 
dead band, °F 

actual pump exit temperature, 170° F 

temperature of the mixture entering the primary regen- 
erator, if it is assumed that there is no heat loss in the 
line connecting the mixing tee and the regenerator, °F 

stack inlet temperature, °F 

steady-state temperature, °F or °R 

time, hr 
initial time, hr 

3 

volume, ft 

3 

volumetric flow rate, ft /min 
3 

specific volume, ft /lb 

work done by an actual cell, Btu or W-hr 

the smaller of X ^mC^ and X ^mC ^ 

cs ^ hs 
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A 2> Ag » r 
y 

AHF 

P 

ahf t 



Subscripts: 

avg 

c/e 

cons 



collection of terms used in the analysis for simplification 

lb m H 2° 

specific humidity, — . H 

m y 2 

pressure correction at constant temperature for heat of 

formation of water, Btu/lb -H_ 

’ m2 

temperature correction at constant pressure for heat of 

formation of water, Btu/lb -H~ 

m 2 

collection of terms used in the analysis for simplification 

efficiency of an actual fuel cell 

efficiency of an ideal fuel cell 

fuel -cell efficiency based on free energy 

fuel -cell efficiency based on heat content 

time constant for fuel-cell transient response, hr 
collection of terms used in the analysis for simplification 

thermal capacity rate of cold-side stream, Btu/hr-°F 
thermal capacity rate of hot-side stream, Btu/hr-°F 
torque, lb^-in. 

fraction of stack exhaust flow rate which bypasses the 
primary regenerator 

average 

condenser exit, primary side 
consumption or consumed 
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cs 

elec 

env 

f 

final 

fw 

gen 

giy 

he 

hs 

htr 

i 

in 

initial 

j 

k 

m 

mix 

p/e 

prod 

pump-sep in 
r 

recir 
recir in 
recir out 

0 


cold side 

2 2 
I R loss in motor, A S2 

environment (fig. 22 and figs. 23(a) to 23(d)) 

force, lbj 

final value 

friction and windage losses, W 
generated heat 

glycol (fig. 22 and figs. 23(a) to 23(d)) 
hysteresis and eddy-current losses, W 
hot side 
heater 

products (fig. 1) 
entering stack 
initial value 
reactants (fig. 1) 

heat transfer quantities of the control volume (fig. 1) 

mass, lb 
’ m 

mixture 

pump -separator exit, primary loop 
produced 

pump- separator assembly input 
pressure- jacket inlet (fig. 13) 
recirculation stream 
recirculation stream entering stack 
recirculation stream leaving stack 


rev 


reversible 


reactant (s) in 
steady state 
fuel-cell stack 

fuel-cell operating temperature, °F or °R (in eqs. (30) 



and (31)) 

theor 

theoretical 

tot 

total 

w 

water vapor 


OPERATING PRINCIPLES AND CHARACTERISTICS 
OF THE FUEL CELL 


rin 

ss 

stack 

T 


The conversion of chemical energy to electrical energy within the fuel cell is 
governed by Faraday's law (ref. 1) which states that during electrolysis, or while a 
voltaic cell is discharging, the passage of 1 faraday (96 500 coulombs) through the 
circuit is accompanied by the oxidation of 1 gram- equivalent weight of matter at one 
electrode and the reduction of 1 gram -equivalent weight of matter at the other elec- 
trode. The fuel cell converts chemical energy into electrical energy without an inter- 
mediate step of conversion into heat. Thus, the irreversibilities associated with the 
conversion of heat into electrical energy are eliminated, and the fuel cell has a higher 
ideal efficiency than the Carnot efficiency of the conventional heat engine. 

From the second law of thermodynamics, the maximum useful work (change in 
free energy) which can be obtained from a chemical reaction at constant temperature 
can be expressed mathematically (ref. 2) as 


AG = AH - T AS 


( 1 ) 


where G is the free energy (Gibbs function), H is the heat content (enthalpy of forma- 
tion), T is the absolute temperature, and S is entropy. The change in the enthalpy 
of formation for the chemical process can be expressed from the heat and mass balance 
(fig. 1) as 


AH =y)m.h. - y^m.h. 

. l l . l ] 
i J J J 


( 2 ) 
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NASA- S- 67- 5544 


m. 

1 

(Products) 



Figure 1. - Fuel-cell heat and mass balance. 


If all the chemical components are at a pressure of 1 atmosphere and at a temperature 
of 77 ° F, then AH = q, where q is the heat of combustion in the conventional combus- 
tion process. However, in the fuel- cell arrangement, q is converted directly into 
electrical energy. The change in entropy of an isothermal chemical reaction is 


AS = Y'm.s. -Vm.s. 

Y 1 1 j J J 


(3) 


The efficiency of a chemical process must be evaluated differently than for the 
conventional heat engine. Efficiency can be defined in two ways. 


_ actual useful work _ power x time 
^AG - maximum useful work ~ AG 


(4) 
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and 


_ actual useful work _ power x time (i 

'AH = change in heat content ~~ AH 

Since AH = AG + T AS, 77 < ?) for the same power output. This relationship is 

An Av_r 

shown in figure 2. 


NASA-S-67-5544 


100 percent 



Maximum 

useful 

work 


Figure 2. - Energy diagram for the fuel cell. 


If the total energy based on the higher heating value (reference 14. 7 psia, 77° F) 
could be converted into electrical energy, then a theoretical potential of 1. 48 V/cell 
(ref. 3) could be achieved, as shown in figure 3. The theoretical potential based on 
the lower heating value for the same conditions is also shown for reference purposes. 
Because of the T AS limitation, the maximum theoretical potential of the cell at 
14. 7 psia and 77° F is 1. 23 volts. This is the voltage which the cell would produce if 
the free energy could be converted entirely to electrical energy with no losses. The 
actual work in an actual fuel cell is less than the maximum useful work because of other 
irreversible phenomena in the physical process. Three types of irreversibilities 
(ref. 3) are encountered which are classified as (1) activation polarization, (2) ohmic 
polarization, and (3) concentration polarization. Not much is known about activation 
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NASA-S-67-5537 



Figure 3. - Hydrogen-oxygen fuel-cell 
performance based on equilibrium 
conditions. 


polarization except that it is attributed to 
various irreversible effects, perhaps 
side reactions at the fuel- cell electrodes. 
Ohmic polarization is caused by the 
electrical resistance of the electrolyte to 
the flow of electrons. Concentration po- 
larization is attributed to concentration 
gradients within the electrolyte (caused 
by an increased water production rate 
at the electrodes) and a possible density- 
gradient effect within the reactant gas 
cavities (attributed to the presence of 
impurities). Because of these losses, 
the performance curve of an actual fuel 
cell would be similar to performance 
curve A of figure 3 for the reaction tak- 
ing place at 14. 7 psia and 77 ° F, with 
reactants and products entering and 
leaving under the same conditions. In 
figure 3, performance is expressed in 
terms of voltage as a function of current 
density, where current density is defined 
as the current output of the cell per unit 
of active cell area. The actual output 
power would be 


P = El (6) 


At higher pressures and temperatures, the maximum theoretical potential in- 
creases because of increases in free energy at higher pressure and temperature. Since 
the Apollo fuel cell operates at a higher temperature and pressure than the standard 
conditions of performance curve A, the performance of the fuel cell will be higher than 
indicated in performance curve A. Although it is a variable -temperature fuel cell, rep- 
resentative single-cell performance is given by performance curve B for 60 psia and 
400° F. 

The efficiency of an ideal fuel cell based on heat content AH (ref. 2) is 


- AG 
^ideal AH 


(7) 


or, by using equation (1) 


^ideal * 


T AS 
AH 


(8) 
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Thus, even an ideal fuel cell operating reversibly and isothermally will have an effi- 
ciency of less than unity, and the heat quantity T AS is exchanged with the sur- 
roundings. 

If all the free energy could be converted to electrical energy in the cell, then the 
maximum useful work would be given by 


AG = E It 
rev 


( 9 ) 


where E rev is the reversible electromotive force of the cell, I is the current output, 
and t is the time required to consume 1 mole of fuel. Ideal efficiency then becomes 


^ ideal 


E It 
rev 

AH 


( 10 ) 


Because of the losses previously described, the waste heat produced in the actual 
cell Q A is greater than the T AS heat of the ideal cell. This is expressed as 


Q a = T AS + Q-^ 


( 11 ) 


where is that portion of the waste heat generated in the actual cell because of 
losses. Thus, the work done by the actual cell is 


W = E It = AH - Q. 

ci A 


( 12 ) 


where E is actual cell voltage. From equation (12), the actual cell efficiency is 
a. 


_ V* _«A 

^actual AH AH 


( 13 ) 


where Q A > T AS. 


The overall chemical reaction of the fuel cell is 


2Hg + Og — 2^0 + heat + electrical energy 


( 14 ) 
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Figure 4 illustrates the chemical balance within the cell which corresponds to this 
overall reaction. At the anode (hydrogen electrode), two hydrogen molecules are ion- 
ized to form four hydrogen ions and four electrons. The electrons travel through the 
external circuit and eventually reach the cathode (oxygen electrode) where they combine 
with one molecule of oxygen and two molecules of water (taken from the electrolyte) to 
form four hydroxyl ions. The four hydroxyl ions then combine with the four hydrogen 
ions from the anode to form four molecules of water at the hydrogen electrode. Two of 
the water molecules go to replenish the potassium hydroxide and water (KOH • HgO) 

electrolyte solution; the other two water molecules leave the cell as steam, having 
absorbed both sensible and latent heat. 


NASA-S-67-5538 


Reactant supply 

I 

I 

\ 

2H„ 


Hydrogen 

electrode 


I 

I 

I 

i 

. C >2 + 2H 2 0 + 4e~ 

t 

I 

I 


Electrolyte 


Reactant 

supply 


From 

external 

circuit 


Oxygen 

electrode 


2H 2 0 


4H 

1 

+ 

4 

I 

l 

1 

40H 


+ 


To external circuit 

t 

I 

1 

4e~ 


2H 2 0, 


Evaporated 

product 


2H 2 + 0 2 


Overall reaction 


2H 2 0 + electrical energy 


Figure 4. - Fuel-cell chemical balance. 


Although the evaporated product (water) takes with it a considerable amount of 
heat on leaving the cell, the waste heat generated by the inefficiency of the reaction is 
greater than that removed by the water vapor. The excess waste heat is removed by a 
recirculating hydrogen stream. The recirculating hydrogen also provides a mechanism 
for water removal. To illustrate better the reactions and their mechanisms within the 
fuel cell, a cell diagram (fig. 5) is presented. 
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Hydrogen enters the hydrogen gas 
cavity where it reacts at the hydrogen 
electrode according to the equation 


H 2 ^ 2H + + 2e~ (15) 


The hydrogen electrode potential depends 
on hydrogen gas pressure and hydrogen 
ion concentration. 

The oxygen electrode reaction is 


0 2 + 2H 2 0 + 4e" C 40H“ (16) 


The oxygen electrode potential depends on 
oxygen gas pressure and hydroxyl ion 
concentration. The total cell voltage is 
the difference between the hydrogen and 
oxygen electrode potentials. 

Two reaction mechanisms exist at the hydrogen electrode. In the first mecha- 
nism, the hydrogen gas adsorbs on the dry nickel electrode surface. The adsorbed 
hydrogen migrates along the nickel surface until it is ionized, whereupon the hydrogen 
ions go into solution in the electrolyte. In the second mechanism, the hydrogen gas 
dissolves in the electrolyte. The dissolved hydrogen diffuses to the nickel surface, 
where it is ionized; the hydrogen ions then go into solution in the electrolyte. The 
second reaction mechanism apparently dominates. 

At the oxygen electrode, the same two physical mechanisms occur as those de- 
scribed for the hydrogen electrode, with the second mechanism again dominating. The 
chemical mechanism is controversial, however, and two theories are advanced. The 
first (the four-electron reaction) is given by equation (16). This theory is favored for 
high-temperature ranges. The second theory (the Berl mechanism) holds that hydro- 
gen peroxide ions are produced at the oxygen electrode and subsequently are decom- 
posed at high temperatures. The chemical reactions corresponding to this theory are 



51.5 psia 

Figure 5. - Fuel-cell electrochemical flow 
schematic. 


20 2 + 2H 2 0 + 4e" - 20H" + 2H0 2 


(17) 


and 


2HO 


catalyst 

2 high temperature* 


20H' + O f 


(18) 
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or, expressed as an overall reaction 


0 2 + 2H 2 0 + 4e" ^ 40H - (19) 


The overall reaction is the same as the reaction in equation (16). The presence of 
hydrogen peroxide ions lowers the electrode potential, thereby contributing to cell 
inefficiency. However, these ions are decomposed more rapidly at higher tempera- 
tures. 

Figure 6 shows schematically the reaction mechanisms occurring at the hydrogen 
and oxygen electrodes. The location at which the electrode, the electrolyte, and the 
reactant gas are all in contact is called the triple interface point. 

NASA-S-67-5541 



(a) Hydrogen electrode. 
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(b) Oxygen electrode. 

Figure 6. - Reaction mechanisms at hydrogen and oxygen electrodes. 



As shown in figure 6(a), the hydrogen gas flows down the coarse pore where (by 
one or both of the two mechanisms previously discussed) it is ionized at the triple 
interface. The hydrogen ions then go into solution in the electrolyte where they react 
with the hydroxyl ions coming from the oxygen electrode to form water. One-half of 
the water evaporates into the hydrogen gas. The hydrogen and water -vapor mixture 
then diffuses out of the coarse pore and is carried away by the hydrogen recirculation 
stream. The remaining water diffuses through the fine pore and into the electrolyte. 
Hydroxyl ions diffuse from the electrolyte, through the fine pore, to the reaction zone 
to replace those used up in the reaction. 

As shown in figure 6(b), oxygen flows through the coarse pore to the triple inter- 
face where it reacts with water from the electrolyte and removes electrons from the 
nickel to become hydroxyl ions. These hydroxyl ions diffuse through the fine pore into 
the electrolyte to replace those used at the hydrogen electrode. Water from the elec- 
trolyte diffuses through the fine pore to the reaction site. 


DESCRIPTION OF FUEL-CELL SYSTEM SELECTED FOR ANALYSIS 


NASA-S-67-5341 


The Apollo fuel cell produces direct- 
current electrical power over a normal 
range of 563 to 1420 watts at a normal 
voltage range of 27 to 31 volts. The mod- 
ule (fig. 7) is 44 inches high by 22. 5 inches 
in diameter and weighs approximately 
245 pounds. Three of these modules, or 
power plants, connected electrically in par- 
allel will be used in the Apollo spacecraft 
to provide electrical power and potable 
water. The module is composed of four 
distinct sections or systems, (1) an energy 
conversion section, (2) a reactant- control 
system, (3) a thermal- control and water- 
removal system, and (4) the necessary 
instrumentation. The last three systems 
are included in the accessory section. 

The energy conversion section is 
shown in figure 8. It consists mainly of a 
stack composed of 31 Bacon-type, series- 
connected cells with associated gas mani- 
folds and connecting leads. The energy 
conversion section is housed in a pressur- 
ized jacket which rests in an insulated sup- 
port assembly. The primary bypass valve 

shown in the figure will be discussed in another section of this report. 



Figure 7. - Apollo fuel-cell module. 


The components forming the accessory section are mounted on a Y-frame as 
illustrated in figure 9. The three legs of the Y-frame are 120° apart. The accessory 
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NASA-S-67-5344 



Figure 8. - Energy conversion section. Figure 9. - Accessory section — Bay I. 


section consists of a nitrogen pressurization system, three regulators, a primary loop 
(hydrogen and water vapor), and a secondary loop (glycol and water), as well as heat ex- 
changers, motor -driven pumps, and plumbing. A condenser connects the two fluid loops. 
Figures 9, 10, and 11 show the various components of the accessory section mounted 
on the three legs of the Y-frame. 

Before examining the system diagram, a discussion of single- cell operation 
would be advantageous. In figure 5, the relative pressure differentials across the 
electrodes are shown. The KOH • HgO electrolyte solution is pressurized by a nitro- 
gen blanket and regulated to 51. 5 + 0. 5 psia. The reactant regulators, with the nitro- 
gen pressure used as a reference, maintain differential pressures of 8. 5 ± 0. 5 psi and 
10. 5 + 0. 5 psi for the hydrogen and oxygen, respectively, above the nitrogen pressure. 
Two parameters governing the performance of the fuel-cell system are the operating 
pressure of the system and the relative pressure differentials across the electrodes. 
The pressure differential across an electrode determines the location of the triple 
interface previously discussed (fig. 6). By extensive testing, the combination of 
pressure and pressure differentials shown in figure 5 has been found to be optimum for 
this system from the combined standpoints of performance and operational feasibility. 


4 
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Figure 10. - Accessory section — Bay II. 
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Figure 11. 


Accessory section — Bay in. 
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Figure 12 illustrates the construction of an actual single cell. The two electrodes 
within each cell are composed of dual -porosity sintered nickel which is formed from 
nickel powder pressed into sheets. The coarse pores (approximately 40 p) are on the 
gas side of the electrode, and the fine pores (approximately 10 p) are on the electrolyte 
side. The two electrodes are similar in construction, but the oxygen electrode has a 
coating of black, lithium -impregnated nickel oxide on the electrolyte side to inhibit 
oxidation. The electrode materials serve as a catalyst in the electrochemical reaction 
and are resistant to corrosion by the electrolyte. A pure-nickel backup plate is used to 
support each electrode and also acts as a gas housing. A Teflon seal, which extends 
around the periphery of the cell, contains the electrolyte and acts as an electrical 
insulator. Although the electrodes are only about 8-1/2 inches in diameter, the entire 
cell is approximately 11-1/4 inches in diameter. The diaphragm section (between the 
electrodes and the cell spacer) accommodates changes in electrolyte concentration as 
the flexible backup plates expand and contract. The 31 cells are stacked in series and 
held together by torsion tie rods. 
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Figure 12. - Single -cell diagram illustrating cell construction. 
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Figure 13 is a schematic diagram of 
the system. Certain components not essen- 
tial to the objectives of this analysis are 
omitted. The diagram is coded to aid in 
distinguishing the different fluid paths. 

The nitrogen subsystem (not shown 
in the diagram) is composed of a small 
nitrogen tank (which holds approximately 
0. 5 lb of nitrogen at 1500 psia), a nitrogen 

regulator, and connecting lines. The regu- 
lated nitrogen pressure (51. 5 ± 0. 5 psia) 
serves a threefold purpose. It is used as 
a reference pressure for the hydrogen and 
oxygen regulators; it is used as a head 
pressure in the glycol accumulator; and it 
pressurizes the jacket around the stack, 
thus pressurizing the electrolyte in each of 
the 31 single cells. 
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Figure 13. - Schematic diagram of fuel- 
system. 


cell 


Hydrogen and oxygen are supplied to 
the module from a cryogenic storage sys- 
tem. The hydrogen is stored at a nominal 
245 psia, and the oxygen is stored at a nom- 
inal 900 psia. The gases are warmed by 
flowing through the connecting lines be- 
tween the cryogenic storage system and the 

fuel-cell system. Then, the gases enter the reactant preheaters before being regulated 
to normal operating pressures. The hydrogen and oxygen subsystems are both equipped 
with purge valves which, when electrically energized, permit a continuous flow of addi- 
tional reactant through the cells. The surplus is dumped overboard. The purging 
process is performed at regular intervals to remove impurities carried into the cells 
by the reactants. 


The makeup (or consumption) hydrogen enters the primary loop at the pump- 
separator exit. There it mixes with the recirculating hydrogen and water vapor and 
proceeds into the pressure jacket through the primary regenerator, where the mixture 
is heated, and from there into the stack. The primary (or hydrogen) loop consists of 
the primary regenerator and bypass control, the hydrogen pump-separator-motor as- 
sembly, a condenser, and an in-line heater for temperature control under low-power 
conditions. 


The primary-bypass-valve sensor detects stack exhaust temperature, which is 
essentially equal to stack temperature. The sensor is a bimetallic strip that also acts 
as a flow diverter. Under high-power conditions when a large amount of heat must 
be rejected, the stack temperature is high, and the bypass valve is open (this is a 
proportional- control valve). Under low-power conditions when heat must be con- 
served, the bypass valve is closed which permits maximum regeneration. 
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The pump- separator is a positive-displacement (vane-type) unit. It circulates the 
hydrogen and water-vapor mixture through the cells to remove waste heat and product 
water. Liquid water from the condenser is separated from the gas stream by centrifu- 
gal action. Input power to the motor (approximately 85 watts) is supplied by three- 
phase, 400-cycle, 115-volt spacecraft inverters. 

The condenser serves a twofold purpose. First, it maintains the primary-loop 
heat balance by rejecting waste heat to the glycol loop for transfer to the radiators. 
Second, it maintains the mass balance in the primary loop by condensing the product 
water vapor from the cells before this water is removed by the separator. 

The secondary loop uses a coolant mixture of ethylene glycol and water. The loop 
consists of a glycol pump, the condenser and preheaters previously discussed, a 
coolant accumulator, and a secondary regenerator and bypass valve. The positive- 
displacement (vane-type) glycol pump circulates the coolant through the secondary- loop 
components and the radiator system. The desired flow characteristics are obtained by 
using a bypass line with an internal orifice to decrease pump flow with increasing 
pressure differential across the pump. Power for the pump (approximately 25 watts) is 
provided by the same spacecraft inverters that supply the hydrogen pump. The coolant 
accumulator maintains a constant pressure within the coolant system regardless of vol- 
umetric changes caused by coolant temperature variations. 

The secondary regenerator controls the heat transferred from the module to the 
spacecraft heat-rejection system to provide the condenser with a relatively constant 
coolant inlet temperature. The bypass valve, which is controlled by the condenser exit 
temperature on the primary side, modulates the glycol flow through the cold side of the 
secondary regenerator. If the primary- side condenser exit temperature becomes too 
high, more of the glycol flow bypasses the secondary regenerator. Less of the glycol 
flow bypasses the secondary regenerator as the temperature goes down. 


GENERAL ANALYTICAL APPROACH 


The basic approach to an analysis of this type involves the formation of an analyt- 
ical model for the fuel cell in a transient mode of operation, with the provision that as 
time t approaches infinity the model must represent the steady -state operation of the 
fuel -cell system. Another stipulation is that at some initial time t^ an initial fuel- 

cell temperature T^ must be known. This initial temperature is an input to the sys- 
tem. As indicated in figure 14, a control volume is set up around the fuel-cell stack. 

By analyzing the incoming and outgoing fluid streams and also the energy released with- 
in the system and transferred from it, a heat balance and a mass balance are obtained 
which form the basis of the analysis. 

Because of the complexity of the system and the interactions of its many compo- 
nents, a completely theoretical approach to the transient problem would lead to a large 
and unwieldy set of simultaneous, higher order, nonlinear partial-differential equa- 
tions. The numerical methods employed for the solution of such a set of equations would 
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almost certainly have stability and convergence problems which are beyond the scope 
of this presentation. 

To circumvent these problems, an analytical approach is undertaken which uses 
power plant and component performance data from various tests conducted throughout 
the development program. This semiempirical approach, used in conjunction with a 
set of simplifying assumptions, aids in holding the number of varying parameters to a 
minimum. Thus, the governing equation for the variation of fuel-cell operating temper 
ature with time has an exact solution. A discussion of the necessary simplifying 
assumptions used in the analysis follows. 


Type of Transient Analysis 

There are two types of transients associated with the operation of this fuel-cell 
system, (1) extremely short-term transients (on the order of milliseconds) and (2) long 
term transients (minutes or hours). The short-term electrical transients, caused by 
the application of "pulse loads" or "spikes, " have come and gone before the fuel-cell 
temperature has had time to change and, therefore, do not contribute to transient 
thermal effects. The long-term transients, however, are associated with a change in 
stack temperature and thus are the type with which this analysis will be concerned. 

To illustrate the long-term thermal transient effect, consider the operating 
characteristics of the power plant on the polarization (voltage -current) curve in 
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figure 1 5. Suppose that the power plant is 
initially at a steady- state operating con- 
dition at temperature Tp current Ip 

and voltage Ep At time t^, a step load 

is applied, so that the current demand 
rises instantaneously from Ij to Ig. 

The voltage at this instant drops from 

E„ to E„ . ... , along the constant tem- 
1 2, initial & 

perature line Tp neglecting the short- 
term electrical transient. The stack is 
now producing a higher current. Hence, 

2 

I R losses are greater, and the stack 
temperature begins to rise at a constant 
current I If this current is maintained, 

temperature and voltage will approach 
steady-state values as t approaches 
infinity. As stack temperature rises, so 
does performance (voltage). This per- 
formance increase with time is illustrated 
in figure 16. Figure 17 illustrates the 
performance variation for a step load 
decrease in current. 
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Figure 15.- Transient fuel-cell 
performance. 
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Figure 16. - Transient fuel-cell perform- 
ance (step current increase). 
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Figure 17. - Transient fuel-cell perform 
ance (step current decrease). 
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Current Variations 


Because of the nature of most spacecraft loads, only step current variations in 
load are considered, with the fuel-cell output current to the spacecraft load being con- 
stant between each step load change. Actually, the current changes slightly because 
of fuel-cell output-voltage variations and because the voltage impressed upon a purely 
resistive load determines the current through it in accordance with Ohm’ s law. 


Fuel-Cell Environment Temperature 

Variations in the environment temperature of the fuel-cell module cause corre- 
sponding changes in the amount of heat lost by the module to its surroundings. These 
changes in heat loss affect module temperature. In this analysis, the fuel-cell envi- 
ronment temperature was a constant 77° F. 


Thermal Mass 

Since the mass of each of the primary- and secondary -loop components is almost 
negligible compared to the mass of the stack, it will be assumed that the thermal 
masses of the components and connecting lines can be neglected. This is equivalent to 
stating that there is no thermal time lag associated with any components other than the 
stack. 


Chemical Equilibrium 

Chemical equilibrium is defined (ref. 4) as that condition in a reversible reaction 
in which the speeds of the reactions in both directions are equal. The chemical equi- 
librium point within the cell depends on many factors, for example, the position of the 
triple interface at the electrodes. Three of the most important of these factors are 
temperature, pressure, and concentration. The analysis accounts for variations of 
temperature, concentration, and partial pressure of the water above the KOH in the 
hydrogen gas cavity. The partial pressure of the water above the electrolyte can never 
reach a state of equilibrium in a transient mode of operation since the concentration of 
the electrolyte depends on partial pressure and temperature, which are always chang- 
ing. The condition of equilibrium is assumed, however, for purposes of calculating 
electrolyte concentration. Furthermore, since the flow rates are small and the flow 
paths are tortuous, it is assumed that the fluid stream leaving the stack is at the aver- 
age temperature of the stack. It is also assumed that the partial pressure of the water 
vapor at the stack exit is equal to the average partial pressure of the water vapor in 
each cell. From this average stack temperature and partial pressure, an average 
electrolyte concentration is calculated using the previous equilibrium assumption. In 
other words, the concentration, temperature, and pressure gradients across the 
31-cell stack are neglected. 
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Stoichiometry 


For the chemical reaction which takes place within the cells, the condition of 
stoichiometry is assumed. The only two conditions which would prevent stoichiometry 
are (1) hydrogen diffusion through the electrolyte and (2) hydrogen or oxygen leakage 
from the system. Hydrogen diffusion through the electrolyte has been shown to be neg- 
ligible by numerous experiments on single cells. Hydrogen leak detectors in the ther- 
mal vacuum chamber have failed to give evidence of hydrogen leakage on power plants 
which have been properly connected to their reactant supply systems, and oxygen con- 
tainment problems are even less severe than those for hydrogen. 


Perfect-Gas Law 

The conditions of pressure and temperature around the primary loop are approx- 
imately 60-psia total pressure (with the partial pressure of the water vapor never 
exceeding 10 psia) at temperatures from 150° to 450° F. Under these conditions, the 
perfect- gas law (ref. 5) is satisfactory and is used throughout the analysis for the hy- 
drogen and water-vapor mixture. 


Condenser- Exit-Temperature Variations 


The variation of condenser exit temperature on the primary side is always less 
than ±5° F from its nominal 160° F, except under adverse glycol-coolant inlet - 
temperature conditions. To account for these extreme variations in primary -side 
condenser exit temperature, a profile of module glycol inlet temperature versus time 
would be needed. Such a profile is not available at this time. For this reason, and 
because the small variations in condenser exit temperature which occur during normal 
operation do little to change condenser performance and cause only a very small change 
in the latent heat of the water condensed h, , a constant 160° F will be used for 


primary- side condenser exit temperature. 


■fg> 


Pressure Conditions 

Although certain very small pressure fluctuations exist from time to time in the 
nitrogen blanket that pressurizes the electrolyte (causing corresponding fluctuations of 
pressure within the hydrogen and oxygen gas cavities), a constant nitrogen pressure of 
52 psia is assumed for the analysis with a constant 8. 5-psid and 10. 5-psid pressure 
differential for the hydrogen and oxygen, respectively, referenced to the nitrogen pres- 
sure as in figure 5. It is further assumed that the pressure drop in the components and 
connecting lines around the primary loop (a total of less than 1 psid) can be neglected. 

A constant total pressure of 60 psia is used for all locations in the primary loop. 
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Primary- Loop Volumetric Flow Rate And 
Efficiency of Water Separation 

For primary-loop calculations, it is assumed that the volumetric flow rate of the 

3 

hydrogen pump- separator assembly is constant at 3. 5 ft /min. Although changing 
condensate mass flow rates at the condenser exit (caused by changing water -production 
rates and variations in condenser performance) cause small fluctuations in torque at 
the motor shaft, bench tests have shown that the pump volumetric flow rate is almost 
constant. These tests have also shown that the efficiency of separation of liquid water 
from the recirculating hydrogen and water-vapor stream is nearly 100 percent. 


Performance Degradation 

Because of the buildup within the cells of impurities that are carried in by the 
reactant gases, the power plant will experience a small degradation of voltage with 
time, under given conditions of current and temperature. Over the operating lifetime 
for which the Apollo fuel cell is qualified (400 hours), this performance degradation is 
statistically less than 1 percent and, therefore, will be neglected in the analysis. 


Reactant Inlet Temperature 

A constant reactant (hydrogen and oxygen) inlet temperature of 77° F will be used 
for the purposes of this analysis, since no profiles of reactant temperature versus 
time are as yet available. It will be evident from the analysis that unless variations in 
reactant inlet temperature become very extreme, performance variations are almost 
undetectable. 


Homogeneity of Hydrogen and Water-Vapor Mixture 

It will be assumed in the analysis that the recirculating hydrogen and water-vapor 
mixture is homogeneous throughout. Although the flow rates are fairly low, the flow 
regime is always turbulent because of the many bends and turns in the primary loop and 
because of such components as the heat exchangers, the in-line heater, and the in-line 
condenser-exit-temperature sensor. Therefore, a homogeneous solution is a realistic 
assumption. At the condenser exit, where two-phase flow exists, the gas mixture 
stream and the condensate stream are treated separately. In the analysis, the forego- 
ing analytical approach is used in combination with the preceding simplifying assump- 
tions and empirical data (where necessary). 


DEVELOPMENT OF GOVERNING EQUATIONS 


The analytical description of the physical behavior of fuel- cell operation is based 
on the proper application of the principles of conservation of mass and energy. Since 
the processes within the fuel cell involve chemical and phase changes and the flow of 
various fluids, the equations must be developed so that the number of independent vari- 
ables is decreased. 
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For developing the reactant-consumption and water-production equations, the 
overall chemical reaction for the formation of water can be written as 


H 2 + 2 °2 H 2° 


(20) 


The mass balance is 


( 2 - 0161b m) H 2 + ( 161b m)°2- 


(18.016 lb m )H 2 G 


( 21 ) 


or 


0 lb m) H 2 + ( 


7. 94 lb 


J m )°2 * ( 8 


94 lb 1H,0 
my 2 


(22) 


With the assumption of stoichiometry, the oxygen consumption and water -production 
rates can be written in terms of the consumption rate of hydrogen as 


rfi~ 

Cleons 


7.94m 


H 2 cons 


(23) 


and 


m n ~ , = 8. 94m TT 

HgO prod HgCons 


(24) 


The term reactant "consumption" is used to distinguish this quantity from total reactant 
usage, which is the sum of reactant consumption and purge (reactants which are vented 
to rid the cells of accumulated impurities). 

On a time basis, Faraday's law can be used to define the electrochemical con- 
stant Kj as 



96 500 A Q( - A- sec 

1.008 gm/sec-H 2 cons gm-H 2 cons 


453. 6 gm/lb 
m 

3600 sec/hr 


= 12 060 A-hr/lb -H 0 cons 
m 2 


(25) 
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Hence 


= 8. 292 x 10 5 lb -H 0 cons/ A-hr 

K. m2 


(26) 


Therefore, by assuming stoichiometry, the hydrogen consumption rate for each of the 
31 cells is 


m , „ 1= 8.292 X 10 5 I 

HgCons/cell 


(27) 


For a 31-cell module, the consumption rate is 


m„ = 2. 571 x 10 3 I 
HgCons 


(28) 


Figure 18(a) is a flow diagram of the 
reactants, products, and recirculating 
stream entering and leaving the fuel- cell 
stack. The fluids of the recirculating 
stream do not enter into the chemical reac- 
tion. The consumption hydrogen enters the 
stack with the recirculating stream, which 
it joins outside the pressure jacket after 
passing through the hydrogen preheater. 
Consumption oxygen enters the stack after 
being conditioned in the oxygen preheater. 


The stack heat loss Q 


stack’ 


electri- 


cal power output P, heat stored Q , waste 

s 

heat generated Q, and the heat entering 
and leaving the stack by means of the recir- 
culating stream are shown in figure 18(b). 

If the gases entering and leaving the stack 
are treated separately, the first law of 
thermodynamics gives 
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Figure 18. - Fluid-flow diagram and heat 
balance for the fuel- cell stack. 


(rfih)jj^ cong + (“ 1 h)Q^ cons ^ + (rc^Hgrecir in + ^^h^O recir in + (^HgCons 


x HF^ 


= (rfih) 


+ (mh) T 


+ (mh) T 


+ 3.413P + Q 


H 2 0 prod out T yii “ v H 2 recir out T V ““ V H 2 Q recir out T T '’‘stack 


(29) 
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where, for convenience, all enthalpies are referenced to 77° F. Kinetic- and potential- 
energy changes across the stack are negligible. The heat of formation of water HF 
will be taken as the lower heating value of hydrogen (since superheated steam leaves 
the stack) and will also be referenced to 77° F. At 1 atmosphere and 77° F (ref. 5), 
HF 77 o= 51 571 Btu/lb m -H 2 . 

Since the chemical reaction within the fuel cell occurs at temperatures near 
400° F and at a pressure of approximately 4 atmospheres, pressure and temperature 
corrections should be investigated. The correction for pressure (ref. 6) is 


AHF = 


p -/ 


m 


(30) 


The correction for temperature (ref. 7) can be written as 


AHF r 


/ 2 

-J 4C p 


dT 


(31) 


where for the reaction, H 2 (gas) + ^ C> 2 (gas) — HgCKgas) 


AC P = C P^ H 2° vapor ) " [pp( H 2 gas ^ + \ C p^°2 gas ^ 


(32) 


Since the reactants and products are treated as perfect gases, = 0. Hence, 

AHF = 0. The correction for temperature AHF T is included in the energy balance. 


The sole purpose of the recirculating stream of water vapor and hydrogen (fig. 18) 
is to remove heat and water. Hence, the quantity Q recir is defined as the heat 

removed by the recirculating stream and can be written as 


v recir 


"(m: 


I (ri^H recir out ' ^^wrecir 


in] + I^^HgO recir out ^^HgO recir in] 


(33) 
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The total energy available from the heat of formation of water (before the temperature 
correction for HF) is 


Q = m TT x HF__o 

gen HgCons 77 


(34) 


By substituting equations (33) and (34) into equation (29), the energy equation becomes 


(,ilh) H 2 cons + (lflh >0 2 cons + Q gen = (lflh) H 2 0 prod + Q recir + 3 ' 413P + Q stack 


(35) 


Introducing the net heat produced in the stack as a result of fuel-cell inefficiency Q 
gives 


Q = Q + (mh)^ + (mh)„ - (rhh) ^ , 

gen HgCons OgCons HgO P rot * 


- 3.413P 


(36) 


If equation (36) is substituted into equation (35) 


Q = Q , , + Q 

stack recir 


(37) 


That is, the net heat generated in the stack is the sum of the stack heat loss and the 
heat removed by the recirculating stream. If equations (23), (24), (28), and (34) are 
substituted into equation (36) 


Q 


= 2. 57 X IQ ' 3 | J [C p(Tin - T o) ] h cons ♦ 7. 94 [5^ - T„] q 
N 3 . 


.cons 


8 - 94 (Pp( T stack- T o)] HO 


+ 51 571) - 3. 413P 


prod 


(38a) 


where the reference temperature Tq is 77° F. From figure 19 and. reference 8, 

the C p av g of hydrogen between 77° and 450° F is 3.44 Btu/lb^- °R. For oxygen 

under these same conditions C = 0.221 Btu/lb -°R; and, similarly, for 

p, avg m ’ ’ 3 ’ 
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Figure 19. - Specific heats at constant pressure for oxygen, hydrogen, 

and water vapor. 


C p, avg ~ °‘ 454 Btu /lb m -°R. Hence, equation (38a) can be written as 


= 2. 57 x 10" 3 I 


3.44(T. n - 77° F) 


H 2 cons 


+ 7. 94(0. 221)(T. n - 77° F) 


Cleons 


8. 94(0. 454)(T„ i „„,_ - 77 


stack ' ‘ 


+ 51 571 


HgO prod 


- 3.413EI 


(38b; 
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by T 


equation (38b) can 


By assuming 


T = T 

Hgin Ogin 


= 77° F and by denoting 


be reduced to 


T 

stack 


Q = 133. 31 - 0. 01043IT - 3. 413EI (39) 

After a step load change (figs. 16 and 17), the fuel-cell current will be constant at I^, 
and equation (39) can be rewritten as 


Q = 133. 3I 2 - 0. 01043I 2 T - 3.413EI 2 


(40) 


To be able to solve the differential equation readily, an ordinary differential equation 
is desired. Therefore, the temperature and concentration gradients across the stack 
are neglected as a first approximation, and the only independent variable is time. 
Equation (40) will be used subsequently in the analysis to obtain a relationship for stack 
temperature as a function of time by expressing Q and E as functions of temperature 
and other variables. 

From figure 18, for a transient thermal condition in the stack 


Q ..+Q = Q+Q . . + Q . . 

recir in ^s recir out stack 


(41) 


where the fuel-cell power output P is included in the expression for Q in equa- 
tion (38a). It is assumed from calculations based on test data that the heat loss from 
the stack varies only linearly with temperature, that is 


Q 


stack 


QA • T + QB 


(42) 


where QA and QB are empirical constants. Now 


Q 


recir out 


Q . . = (rfiC AT) 
recir in \ p / 


Hgrecir 


+ (* c p 


AT) 


HgO recir 


(43) 


where 


AT 


rji ^ fji _ rp _ rp 

stack out stack in ~ stack stack in 


T - T 

stack in 


(44) 
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If a constant condenser exit temperature of 160 F, a constant total pressure of 

3 

60 psia, and a constant volumetric flow rate of 3. 5 ft /min for the positive- 
displacement-vane pump at the condenser exit are assumed, the perfect-gas law may 
be applied to the fluid stream (ref. 9) as follows 


pV = (mR) 


. T 
mix 


(45) 


therefore 


(mR) 


mix 


(60 X 144) X (3. 5 X 60) 
620 


= 2925 BtuAr/deg 


(46) 


But 


^ rflR ^mix = + ^ rflR ^ 


h 2 ° 


(47) 


therefore 


X 767 Ay „ . x 85. 3 

2925 H 2 recir H 2° recir 


m 


Hgrecir 


m 


Hgrecir 


m 


E^recir 


(48) 


or 


2925 


m 


= 767 + 85. 3 y 


F^recir 


c/e 


(49) 


H 2 0/lb 

perfect-gas equation of state for each constituent of the mixture and by using Dalton's 
law of partial pressures 


dry Hg). By using the 


where Y c j e = specific humidity at condenser exit (lb m 


m 


w 


P V 
R T 


w 


(50) 
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and 



V 

V 


(51) 


where the subscript w refers to water vapor and the subscript Hg refers to dry hydro- 
gen. Now 


m 


R 


r = 


w 


w 


H r 


P. 


R/MW 


w 


H r 


p X MW 
*w w 


P„ 


m 


Hr 


P 


Hr 


R 


w P H 0 R/MW 


w 


p„ x MW 
H 2 H 2 


Hr 


2. 016 


= 8. 94 


P 


Hr 


(52) 


But 


P 



+ P- 


w 


(53) 


therefore 



p = 60 - p 


’ W 


w 


(54) 


or 


y= 8. 94 


■ w 


60 - p 


(55) 
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or, for later use, equation (55) can be written as 


p w “ 


60 


1 + 


8. 94 


(56) 


For a condenser exit temperature of 160° F, the partial pressure of the saturated 
water vapor is 4. 741 psia, hence 


y c/e 


8. 94 x 


4.741 
60 - 4. 741 


(57) 


or 


y , = 0. 767 

r c/e 


(58) 


If equation (58) is substituted into equation (50) 


2925 

^H^ecir 


767 + 85. 3 x 0. 767 = 832. 5 


(59) 


or 


m 


Hgrecir 


3. 51 lb /hr 
ml 


(60) 


therefore 


m. 


HgO recir 


= y / x m„ 
'c/e H r 


recir 


= 0. 767 x 3. 51 


(61) 


or 


rh . = 2. 70 lb /hr 

H 2 O recir m/ 


(62) 
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If equations (60) and (62) are substituted into equation (43), for the average values of 
specific heat stated before, and if stack inlet temperature is denoted by T 


O - O 

recir out recir in 


(niC AT\ + (mC AT") 

' ^ 'Hgrecir ' ^ 'HgO 


recir 


= (3. 51)(3.44)(T - T g ) + (2. 70) (0. 454) (T - T g ) 


(63) 


or 


Q , - Q . . = 13. 3 fT - T ) 

recir out recir in s) 


(64) 


The thermal mass of the stack M from test data is 


M = 32 Btu/° F 


(65) 


Hence, the heat quantity stored in the stack by raising the temperature of the stack 
from Tj to T is 


q g = y’ T MdT = M dT = M(T - 

•'T 1 



( 66 ) 


or 


dq s_ 0 _ dT 

~ar - Q s - M "cir 


(67) 


Finally, if equations (42), (64), (65), and (67) are substituted into equation (41) 


dT 


Q = 32 ~ + 13. 3(T - T g ) + (QA • T) + QB 


( 68 ) 
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or 


Q = 32 + (13. 3 + QA)T - 13. 3T g + QB 


(69) 


If equation (67) is substituted into equation (39) 


32 + (13.3 +QA + 0.010431)T - 13. 3T + (QB - 133. 31 + 3. 413EI) = 0 (70) 

cit s 


From appendix A, if the time degradation of performance caused by the buildup of 
impurities is neglected, the variation of module output voltage with temperature, 
current, and electrolyte concentration is 


E = (CAHI • T + CBHI)I + (CDHI • T + CEHI) - 0. 2473 (percent HgO - 27. 0) (71) 


where 360 0 F < T < 480 ° F, 5 amperes < I < 55 amperes, and 22 < percent HgO < 32. 


For the current region from open circuit to 5 amperes, the equation corresponding 
to equation (71) is 

E - (CALO • T + CBLO)I + (CDLO • T + CELO) - 0. 2473 (percent HgO - 27. 0) (72) 

where 360 0 F s T < 440 ° F, 0 < I < 5 amperes, and 22 < percent HgO < 32. Typical 
empirical values of the voltage curve fit constants are given in the following table. 


High- current region 

CAHI = 7. 96 x 10“ 4 V/A- °F 
CBHI = -0. 481 V/A 
CDHI = 0. 026 V/ °F 
CEHI = 23. 0 volts 


Low-current region 

CALO = 2. 775 x 10 -3 V/A- °F 
CBLO = -1. 903 V/A 
CDLO = 0. 0157 V/°F 
CELO = 30. 24 volts 
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Equations (70), (71), and (72) are the principal equations that describe the tran- 
sient behavior of the fuel-cell system. Other expressions will be introduced, as 
needed, in solving these three equations for appropriate boundary and initial conditions. 


SIMPLIFIED TRANSIENT ANALYSIS 


A simplified transient analysis is an analysis of only the fuel-cell stack and ne- 
glects the various heat exchangers and other parts of the system. Such an analysis is 
developed as a first step toward mathematical formulation of the transient fuel-cell 
problem prior to incorporating the auxiliary components into the analysis. Although 
the results of the simplified analysis are not in good agreement with test data (espe- 
cially at higher temperatures), they do provide satisfactory results within limits. 

The energy equation for the system is equation (70). For the simplified analysis, 

however, the temperature T at which the hydrogen and water-vapor stream enters 

s 

the stack (fig. 13) is assumed to be a constant 350° F; therefore, in equation (70) 
dT 

s 

— gj- = 0. This is equivalent to assuming the presence of a heat exchanger between the 

condenser exit and the stack inlet which always raises the temperature of the stream 
flowing through it to 350° F, regardless of other changing conditions. With this as- 
sumption, equation (70) becomes 


32 ^ + (13. 3 + QA + 0. 010431)T + (QB - 4655 - 133. 31 + 3. 413EI) = 0 (73) 


From equation (71), if variations of electrolyte concentration with time are neglected 
and a constant of 27 percent water by weight is assumed for the electrolyte solution, 
the variation of voltage with current and temperature is 


E = [(CAHI)T + CBHl]l + (CDHI)T + CEHI 


(74) 


where 360 ° F ^ T ^ 480 0 F, 5 amperes < I < 55 amperes, and percent HgO = 27. 0. 

Note that the low-current region from open circuit to 5 amperes is excluded in 
the simplified analysis. If equation (74) is substituted into equation (73) with 
QA • T = 165 Btu/hr = Q gtack and QB = 0 


A -gj- + BT + C = 0 (75a) 
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where 


A = 32 


B = (3. 413CAHI)r + (3. 413CDHI + 0. 01043)1 + 13. 3 


C = (3.413CBHI)! + (3.413CEHI - 133.3)1 - 4490 


(75b) 


Equation (75a) is solved by separation of variables and can be evaluated by applying 
the boundary condition that T = at t = t^. The result is 




(76) 


which is of the form 


T = a + be 


-At 


(77) 


where 


a = - 


C 

B 


(78) 


At 

b = ( T i - a)e (79) 

A = | (80) 

The fuel- cell stack temperature T 
as t becomes infinite; that is, as 

into equation (78), a = T . Thus, 

ss 

of a and b into equation (77), the 


must approach its steady-state temperature T gg 

t — T — T . When this condition is substituted 
ss 

At 0 

b = ^Tj - T gs )e . By substituting these values 
equation for the change in stack temperature with 
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time becomes 


or 


T - T ss + ( T i - T sa > X ‘° • e ' xt 

-X(t-t ft ) 

T = T - (T - T,)e 
ss \ ss 1J 


(81) 


(82) 


If the typical voltage curve fit constants previously stated are used, the stack temper- 
ature T from equation (79), where a = T , is given by 


T = 
ss 


1. 6421 + 54. 81 + 4490 


2. 717 X 10" 3 I 2 + 9. 91 X 10 -2 I + 13. 3 


(83) 


and the steady-state voltage is computed as 


E gg = (7. 96 X 10 -4 T gs - 0. 481)1 + (0. 026T gs + 23. 0) 


(84) 


The steady-state power level of the fuel cell is then given by 


P = E xi 

SS SS 


(85) 


The results of this simplified analysis are compared to those of the more general 
transient analysis in the section on comparison of analytical results. 


GENERAL TRANSIENT ANALYSIS 


The analysis is extended to the general situation with the incorporation of the 
primary-bypass-control and primary-regenerator characteristics into the system. 
Although this analysis will be termed a general transient analysis, it is only general 
in comparison with the simplified analysis and is still subject to the simplifying as- 
sumptions in the section on general analytical approach. The incorporation of the by- 
pass control and the primary regenerator produces a feedback effect in the temperature 
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region where regeneration occurs as will be shown later by comparing the results 
of the simplified analysis to the results of this more general analysis. This feedback 
effect quickens the thermal response of the system considerably in this region and 
narrows the operating-temperature range of the system. 

In this analysis, the following assumptions previously used in the simplified anal- 
ysis will be eliminated. 

1. The temperature T at which the hydrogen and water-vapor mixture enters 

s 

the stack, previously assumed to be constant, will be allowed to vary in accordance 
with the response of the primary bypass valve and the primary regenerator. 

2. The previous assumption of constant electrolyte concentration will no longer 
be used. The concentration of the KOH electrolyte solution will be allowed to change 
according to variations in electrolyte (stack) temperature and in partial pressure of the 
water vapor above the KOH at the gas-electrode-electrolyte interface. 

The assumption of pressure and temperature equilibrium at the gas-liquid interface will 
be retained for the purpose of calculating electrolyte concentration. It will further be 
assumed (as with stack temperature) that the partial pressure of the water vapor at the 
stack exit is equal to the average partial pressure of the water above the KOH within 
the stack. 

In addition, the analysis will include the following considerations. 

1. The effects of the heat input to the mixture stream at the hydrogen pump, 
caused by the inefficiency of the pump, are considered. 

2. Thermal effects caused by the cycling of the in-line heater at low-power 
levels are included. 

3. The low-current region (0 to 5 amperes) is included in the analysis to pro- 
vide the capability of analyzing the hot-standby or no-load condition. 

All previous assumptions, such as a constant condenser exit temperature of 160 ° F 
and a constant fuel-cell-system environmental temperature of 77 0 F, will be retained. 

In the schematic diagram (fig. 13), the heat input to the system at the hydrogen 
pump is calculated as follows. The pump- separator motor is a three-phase, 
Y-connected motor that draws approximately 0. 76 ampere of line current, which is 
equivalent to 0. 76 A/phase. The spacecraft inverter supplies the motor with 
115 V/phase; therefore, the motor is supplied with 87. 3 V-A/phase or a total of 
262 V-A for the three phases. At a nominal power factor of 0. 32, the power input to 
the motor is approximately 84 watts. 

Some of this power input is dissipated in the motor because of losses, the 
remainder being shaft output power. This can be expressed as 


H 


• = P T 

.pump m L 


+ P 


shaft 


( 86 ) 
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where = power dissipated in motor because of losses, and = shaft output 

power of motor, or P g haf t = i 
into three general categories, Pj, 


power of motor, or P g haf t = shaft horsepower x 746. The power losses can be divided 


fw’ 


P, , and P . , where 

he’ elec’ 


P T = P* + P, + P i 
L fw he elec 


(87) 


Since the motor torque is 32. 5 oz-in. at 2125 rpm 


shaft horsepower = gg q q q 


where r = torque ^lb^-in.) , and n = motor speed (rpm), or 


( 88 ) 


x 2125 

shaft horsepower = — gg ggg — 

= 0. 0684 (89) 


so that P sllia ft = 51 watts. This means that P^ = 33 watts, or that the motor efficiency 

defined as the ratio of shaft output power to total input power is approximately 61 per- 
cent. The 51 -watt equivalent of shaft horsepower operates the pump and the water sep- 
arator. 

The losses associated with the pump and the water separator can also be calcu- 
lated. The input power to the motor can then be expressed as the sum of all these 
losses, plus the flow work required to pump the mixture stream around the primary 
loop, plus the work required to impart an amount of kinetic energy to the liquid water 
sufficient to remove it from the system by centrifugal force. 

Regardless of the amount of the total input power converted to heat by losses, 
either directly or indirectly, within the pump- separator assembly, the remainder of 
input power will also eventually appear in the system (although not entirely in the pump- 
separator assembly) as heat from friction. This heat from friction, which is dissipated 
in other parts of the system, is only a very small part of the overall power input. It 
will be assumed that the total input power to the pump-separator assembly is dissi- 
pated as heat in that same assembly. Not all of this heat enters the recirculating 
stream, however. The fraction which is absorbed by the stream will now be calculated. 

As indicated in figure 13, the temperature of the mixture stream at the pump exit 
is found from thermal-vacuum test data (taken from in-house performance evaluations 
of Apollo fuel-cell systems) to be approximately 170 ° F when the fuel-cell system is in 
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a nominal 77 ° F environment. If the small preheating effect is neglected, the reactant 
(makeup) hydrogen enters the primary loop at the mixing tee at a temperature of 77° F. 
If the specific-heat values and flow rates of the hydrogen and water vapor from the sec- 
tion on development of governing equations are used, the temperature to which the 
stream would rise if the total power input were absorbed by the stream is calculated 
from the following energy balance. 


S 


•ump-sep m = ^{ AC p ) ( T p/e ‘ 160 ° F ) 

(^Hgrecir * *"p, Hg + recir X C p, H 2 o) ( T p/e " 160 F ) ^ 


But 


% 


ump-sep m 


= 84 x 3. 413 


= 287 Btu/hr 


(91) 


therefore 


n _ 287 

p/e theor (3. 51 X 3.44) +12.70 x 0. 454) 


+ 160 = 181.6 F 


(92) 


Since the actual temperature (from test data) to which the recirculating stream rises 
as it flows through the pump- separator assembly is only 170° F, the heat loss at the 
pump from radiation and conduction is given by 



£(mC p ) (181. 6 - 170. 0) = 13. 296 x 11. 6 = 154. 1 Btu/hr 


(93) 


Hence, the actual heat input to the recirculating stream at the pump is 132. 9 Btu/hr 
or only about 46. 3 percent of the total input power. 

If the notation in figure 13 is used, the mixing equation for the hydrogen makeup 

tee is 


TI X C „ x (T . - T n ) + m-j . x C TT (T , - T/) + m TT „ 
HgCons P> H 2 ' rin Hgrecir P, Hgl p/e 0) HgO recir 

X C p, H 2 o( T p/e “ T o) = m H 2 recir X C p, H 2 ( T R ’ T o) + m H 2 0 recir 

X C p, H 2 o( T R “ T 0) + m H 2 cons X C p, H 2 ( T R " T o) ^ 
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Ill 


where T n is a reference temperature; T_ is the temperature of the mixture entering 

the primary regenerator, if it is assumed that there is no heat loss in the line connect- 
ing the mixing tee and the regenerator; and T is the actual pump exit temperature 

(170° F). All terms containing the reference temperature Tq are removed by sub- 
traction. K equation (94) is rearranged and if the known temperatures, the flow rates, 
equation (28), and the specific heats are substituted into the equation, the expression 
for T 0 is given as 

rv 


T = 
R 


0.6811 + 2261 
8.841 x 10" 3 I + 13.3 


(95) 


where I is fuel-cell current. It is assumed that T R changes instantaneously with a 

step current change, although there is some small timelag associated with the mixing 
tee. 


To calculate the mixture temperature at the cold- side exit of the primary regen- 
erator T , the cold- side effectiveness E of the regenerator must be considered, 
s OS 

The effectiveness of this counter-flow heat exchanger is defined as 


E s 
cs 


E e%) 


cs 


( T s - T r) 


X (T - T r ) 


(96) 


and 


^( mC p) cs ( mC p, H 2 recir) + ( mC p, HgCons) + ( mC p, H 2 0 recir) (97) 

where S/mC \ is the thermal capacity rate of the cold- side stream, X is the 
V P/ cs 

smaller of £/mC ) and Z)/mC ) , and 22(mC) is the thermal capacity rate 

v p 'cs v p/ hs v p/ hs 

of the hot- side stream. 


To determine X, the characteristic of the primary bypass control must be intro- 
duced. Figure 20 is a plot of percent bypass versus temperature. The bypass valve is 
biased to stack temperature T by means of a bimetallic strip which is placed within 
the valve in the stack exhaust stream. Curves A and B in figure 20 illustrate the 
effects of hysteresis inherent in the valve. For example, if the stack temperature is 
at 430° F and rising, then the percent bypass follows curve A. However, if the stack 
temperature starts decreasing after reaching 444° F, the valve will continue to bypass 
50 percent of the total flow around the regenerator until T = 441° F. At this point, the 
percent bypass will start to decrease along curve B. 
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Stack temperature, °F 


Figure 20, - Primary-bypass- valve characteristics. 


Since the valve-temperature dead band caused by hysteresis is relatively small, 
it is neglected in this analysis. Also, since the normal operating-temperature range of 
the power plant is 380° to 440° F from minimum to maximum power, the linear-bypass- 
valve characteristic (typical) of curve C is used. 

Now 


XYmC 


(mC ) 

v p 4s 


= ( m TT . + m rT \ C TT + ( m TT _ . x C _ 

Hgrecir HgConsJ p, Hg ( HgO recir p, HgO 


= (3. 51 + 2. 57 X 10“ 3 l) x 3. 44 + (2. 70)(0. 454) 


(98) 


or 


L(mC\ = 8.841 x 10 _3 I+ 13.3 
' P ' cs 


(99) 
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Let the bypass fraction be denoted by <p; that is, = percent bypass x 10 , so that 


+ • C „ rTl (100) 

( HgO recxr HgOJ p, HgO 


where _ is that portion of the water produced which exits from the stack. This 

h 2 o 

term can be expressed as 


m H 2 0 m H 2 0 prod + V dt 


( 101 ) 


HO 

where — < 0. In a steady-state condition, all the water produced leaves the 
dm H 2 0 

stack and — g^ — = 0. Therefore, for steady-state operation 


m H 2 Q m H 2 0 prod 


In the transient condition, — dtf~~ ~ But O chan & es as electrolyte concentra- 

2 

tion (percent HgO) changes; therefore 


dlfl H 2 0 d (percent HgO) dm H 2 0 
dt dt ~ d (percent HgO 
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However, since changes in electrolyte concentration are very small, this term can be 
neglected with only very small error; hence 


~ prod 


(104) 


The expression for the hot- side thermal capacity rate is then 

X^mC \ = (1 - 0)f7m tJ . x C „ \ + /rfi„ „ . + m TT ,\ x C TT „ 

V pj hs ^ H 2 recir p, H 2 0 recir HgO prodj p, HgO 

= (1 - 0)|(3. 51 X 3.44) +(2.70 + 8.94 x 2.57 x l(f 3 l)(0.454)] 

= (1 - 0)(lO. 431 x 10" 3 l)+ 13. 3 (105) 


or 

T,(mC\ = (1 - 0)a 

^ p 'hs 


(106) 


where a = a a) = 10. 431 x 10 _3 I + 13. 3. 

It is evident from the graph in figure 20 that there are two distinct stack- 
temperature regions to be considered; that is, the regions of variable regeneration 
from 360° to 420° F and from 420° to 440° F. (The upper limit used is 440° F since 
this temperature roughly corresponds to the temperature of the stack at its rated 
maximum-continuous-output power of 1420 watts. ) It is also evident from appendix A 
that there are two current regions. Because of the various regions of current and 
temperature, there are five cases to be considered. These cases are listed in table I. 
The cut-off temperatures, 420° and 440° F, used in table I are typical; however, these 
values are generalized for computer programing purposes by using TgpjjQ anc * 

T BPHF res P ectivel y- 
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TABLE I. - CASES CONSIDERED IN THE DEVELOPMENT OF THE 
TRANSIENT EQUATIONS FOR DIFFERENT CURRENT AND TEMPERATURE 

REGIONS OF OPERATION 



Current region, 
A 

Temperature region 

In- line 
heater 

Case 



Regeneration 





0 to 5 

5 to 55 

Full, 

0=0 

Partial, 

0=0 

Temperature range, 
°F 

On 

Off 

I 


X 

X 


360 < T < 420 


X 

II 


X 


X 

420 < T < 440 


X 

in 


X 

X 


360 < T < 420 

X 


IV 

X 


X 


T. . , < T < 420 
high 


X 

V 

X 



X 

420 < T < 440 


X 


The transient equations for all cases for the variation, after a given step load 
change, of stack temperature with time are derived in the following sections. 


Case I 

The conditions for case I are the following: 

1. 5 amperes s I < 55 amperes 

2. 360° F < T < T BpLO 

3. Full regeneration ( 0 = 0) 

4. In-line heater off 

To determine X for use in equation (96), the values of X^mC^ and 
S/mC 'l are compared from equations (99) and (105), where in equation (105), 0=0. 

V p 'hs 

At I = 5 amperes 


E(mC \ = 13. 344 Btu/hr-°F 

' P 'cs 


(107) 
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and 


£YmC \ = 13.352 Btu/hr-°F 

' P 'hs 


(108) 


At I = 55 amperes 



13. 786 Btu/hr- 0 F 


(109) 


and 


£/rnC ^ = 13. 874 Btu/hr- °F 

V IV hs 


( 110 ) 


Therefore, 0. 008 < 


£ (“ c p) 1is - s (* c p) 


cs 


< 0.088. It is evident from these calcu- 


lations that the difference in hot-side and cold-side thermal capacity rates is always 
less than 0. 1 Btu/^ir-°F in the region of full regeneration, with £YmC \ only 

\ vhs 

slightly greater than £/mC \ . Therefore, it will be assumed that 

\ p / cs 

£/mC \ = S(niC ^ ; thus, the value of X is chosen to be E(mC \ . The equation 

' P 'hs ' IV cs ' P ' cs 

for cold- side regenerator effectiveness is then simplified considerably to 


E 


cs 


T - T 
s R 


R 


( 111 ) 


By using a constant cold-side effectiveness of 0. 84 in equation (111), the inlet 
temperature to the stack is 


T = E T + (1 - 
s cs \ 



( 112 ) 


By substituting equation (111) into equation (64) and denoting Q recir out - Q recir in by 

O • , the heat removed from the stack by the recirculation stream becomes 
recir’ 


Q . = 13. 3(1 
^recir 


E cs)( T - T R) 


(113) 
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To calculate the steady- state temperature of the fuel cell as a function of current, 
the steady- state heat balance is 


Q " Q recir + Q stack (114) 

dT 

where Q is given by equation (69) with = 0. Substituting equations (69), (42), and 
(113) into equation (114) gives 

- T R )J+ (QA)T + B (115) 

where E is given by equation (71) (repeated here for convenience) as 

E = (CAHI • T + CBHI)I + (CDHI • T + CEHI) - 0. 2473 (percent HgO - 27. 0) (71) 

As evidenced in appendix A (eq. (A6)), percent HgO is given by 

percent HgO = 61. 78 - 0. 0997T + p^ (116) 


1(133. 3 - 0. 01043T - 3. 413E) = 13. 3^1 - E^^ 


where p^ is the partial pressure of the water above the KOH, which is assumed here 
to be the equilibrium pressure calculated at the stack exit. From equation (56) 


P 


w 


60 

8. 94 


(56) 


where y is the specific humidity 
By definition 


lb H 0 0 
m 2 

lb dry H„ 
m J 2 


of the mixture stream at the stack exit. 


y = 


m TT _ . + m TT „ , 

HgO recir HgO prod 


m 


Hgrecir 


(117) 
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or 


y = 


2.70 + 22.98 X 10" 3 I 
3. 51 


(118) 


By substituting equation (56) into equation (116) and then substituting equation (116) into 
equation (71), the voltage relation (eq. (71)) becomes 


E = [(CAHI)T + CBHl]l + [(CDHI)T + CEHl] + 0. 02466T + T (119) 

where y is given by equation (73), F = T(I), and 


r(l) = - 8 . 601 - 


14. 838 

l + ^4 


( 120 ) 


By substituting equation (119) into equation (115), simplifying, and then rearranging, 
the steady-state stack temperature T becomes 

So 


T = 
ss 


(-3. 413 • CBHI • I 2 ) + [133. 3 - (3. 413 • CEHl) - (3. 413 • T)]l + [l3. 3(1 - - 

(3. 413 • CAHI • I 2 ) + [ 0. 09459 + (3. 413 • CDHI)] I + Jl3. 3 (1 - E cg jj + QA 


QB 


( 121 ) 


This steady- state temperature will be used later as a boundary condition. The steady- 
state voltage for this case is given by equation (71) with T = T . For the transient 

o o 

condition, by substituting equations (119) and (112) into the governing differential equa- 
tion (eq. (70)), a linear first-order equation of the form of equation (75a) is obtained. 
However, the coefficients for this case are given by 


A = 32 (122) 

B = (3. 413CAHI)I 2 + (0. 09459 + 3. 413CDHI)I + 13. 3 + QA - 13. SE^ (123) 

c s 
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(124) 


C = (3. 413CBHI)I 2 + (3. 413CEHI + 3. 413 T - 133. 3)1 
+ QB-13.3(1-E cs )T r 

The temperature-time relationship for this case is the same as that of equation (82), 
where T gg = - §. 


Case II 

The conditions for case II are as follows: 

1. 5 amperes < I < 55 amperes 


2. T < T < T 

X BPLO 1 BPHI 


3. Partial regneration 

4. In-line heater off 

Since E/mC \ and E(mC \ are approximately equal with full regeneration, it i 
' P 'cs ' "hs 

evident that with any amount of regeneration E/mC \ will always be less than 

V p /hs 

E(mC \ . Therefore, X = XYmC ^ , and equation (96) becomes 

V p /r.s \ p 'hs 


is 


E 


E 


cs 


= 0 . 84 = 


(Ncs • C 


T - - t r) 


2rp)hs ^ 


r) 


( 125 ) 


Now, the linear primary-bypass-valve characteristic (fig. 20) can be expressed 
mathematically as 


percent bypass = (BP1 • T) + BP2 


(126) 


Typical bypass values of BP1 and BP2 are 2. 167 percent/°F and -910. 4 percent/°F, 
respectively, for Tgp^Q = 420° F. 
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By substituting equations (99), (105), and (126) into equation (125), simplifying, 
and solving for T , it is evident that 


where 


T S = T R + f>l T * W * ^ 


(127) 


-E • a 


fficy (bp1 x 10 ' 2) 

\ v es 


(128) 


E * a 

Q — CS 
Pr, - xr^y- 




[(l - BP2 x 10' 2 ) + (bPI X 10" 2 )t^J 


(129) 


0 * = 


" EcS ' a ' T — [l - ( B P2 X 10 -2 )J 


S’ E(mC p ) 


cs 


(130) 


By substituting equation (127) into equation (64) 


Q recir = (" 13 - 3/3 l) T + [ 13 ‘ H 1 ~ ^) T ] “ 13 ’ 3 ( T R + *a) < 131 > 


or 


Q . = 6 T 
recir 1 


+ 6 2 T 


+ &€ 


(132) 


where 5^ = -13. 3/3^, = 13. 3^1 - /3^j, and 6^ = -13. 3^T^ + /3^. Again, the steady- 

state heat balance as given by equation (114) is 


Q - Q 


recir 


+ Q 


stack 


(114) 
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Substituting equation (69), with = 0, and equation (130) into equation (114) and then 

solving for the steady-state fuel-cell temperature T gives 

ss 


T 


ss 



(133a) 


where 


£ 2 = (3. 413 • CAHI)I 2 + (3. 413 • CDHI + 0. 01043)1 

+ + QA ^ 

= (3.413 • CBHI)I 2 + (3.413 • CEHI - 133.3)1 
+ + QB 


(133b) 


An order-of-magnitude analysis for the current region under consideration reveals that 
the positive root must be chosen in equation (133a) for T to make the mathematical 

SS 

model compatible with the actual fuel cell. Therefore 


T 


ss 


“^2 + V ^2 " 4 ^ 1^3 

^1 


(134) 


The order-of-magnitude analysis also insures that for the current region considered 
2 

the discriminant ^ a l wa Y s greater than zero, thereby avoiding the possi- 

bility of imaginary roots. 

By using equation (71) in conjunction with equation (134), the steady-state voltage 

for this case can be obtained. The differential equation for the transient condition is 

derived by substituting equations (119) and (127) for E and T , respectively, into 

s 
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equation (70). This yields 


32 § + 5 1 t2+ (V 6 4) T - C 1 = ° 


(135) 


where 


6 4 = -(3.413* CAHI)I 2 - (3.413 • CDHI + 0. 01043)1 - QA (136) 


and 


C 1 = -(3.413 • CBHI)I 2 - (3.413 • CEHI - 133.3)1 - 6 g - QB (137) 


/dT \ 

Note that for steady- state conditions — Oy , the constant can also be evaluated 
as 


C 1= 6 1 T ss 2+ (V 5 4) T ss < 138 > 

where T is given in equation (134). 
s s 

Equation (135) can be rearranged and integrated as 



32 dT 

' 2 
a + bT + cT 



(139) 


2 

where a = C^, b = (S 4 - 6 2 ), and c = -6^. Now, by defining q as 4ac - b , the 
value of the stack temperature at any time after a step load change from 1^ to I 2 
(where the stack temperature was T^ at t^ when the load changed instantaneously) is 


T = - 




tan 


-1 


m 



(140) 
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for q > 0, and 


(141) 



for q < 0, where the positive root is always used for and V^q. 
The voltage at any temperature T is given by equation (71). 


Case III 


The conditions for case III are as follows: 

1. 5 amperes ^ I < 55 amperes 

2. 360° F < T < 420° F 

3. Full regeneration (<p = 0) 

4. In-line heater on 


Although the equations for this case are good for 360° F ^ T < 420° F, the in- 
line heater is designed to sense stack temperature, so that if the stack temperature is 
falling the heater is activated when the temperature reaches some temperature T, 


in the range of 380° ± 5° F. 

risen to some temperature T^j ^ in the range of 390" + 5 


low 

The heater then stays on until the stack temperature has 

F. The aptual values of 


low 


and 


T^gh var y ^ rom P ower plant to power plant. 


The in-line heater dumps 


heat into the power plant for low-power operation according to the equation 


Q 


htr 


3.413E 2 

^tr 


(142) 


Since the voltage varies only slightly within the relatively narrow temperature dead 
band (only 3° or 4° on most power plants) in which the heater is on, a constant heater 
power based on 31 volts will be used in the analysis. Heater resistance is 6. 02 ohms. 
By substituting these values into equation (142), the heater power is found to be 
159. 6 watts, or 544. 8 Bfu/hr. It is assumed that all of this heat is used in raising the 
temperature of the power plant. For the purpose of calculating a steady- state temper- 
ature with the heater on, equation (36) becomes 


Q 


+ ^htr ^ stack 


+ Q 


recir 


(143) 
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By using the voltage equation (eq. (119)) in conjunction with (1) equation (140) for Q 
dT 

(with = 0), (2) equation (113) for Q recir , (3) Q stack = (QA • T) + QB, and 

(4) Q^ r = 544. 8 Btu/hr, the steady- state stack temperature with the heater on is given 

by 


13. 3(1 - E cg )T R + 544. 8 - QB + 133. 3I tot |^CBHI • ^ + CEHI + r“| 

T ss = 1373(1 - E cs y + 0. 01043I tot + 3. 413I tot |pAHI "• 1^ +^CDHI + 0. 02466 j +T3A 

(144) 

where T is calculated on the basis of total current I, ,, and 1=1 + 1, =1 + 31/6. 02, 

tot tot htr ' ’ 

or 


I, . = I + 5. 15 
tot 


(145) 


The transient equation for this case is the same as equation (82) in case I, except 

that the steady- state temperature T is the same as in equation (144). 

s s 


Case IV 


The conditions for case IV are as follows: 


1 . 

2 . 


0^is5 amperes 


T •< T < T 
high - 1 - a BPLO 


3. Full regeneration (<p = 0) 

4. In-line heater off 


To calculate T for case IV, equation (72) is modified by substituting equation (56) 
ss 

into equation (116) and then substituting equation (116) into equation (72) to obtain 


E = [(CALO • T) + CBLOjl + (CDLO + 0. 0247)T + T 


(146) 


where r =-8.601 


14. 838 

77E5' 
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For case IV, the steady-state temperature, obtained by substituting equation (146) 
into equation (115), is 


T = 
ss 


(-3. 413CBLO)! 2 + (133. 3-3. 413CELO - 3. 413T)I + 13. 3(1 - E cg )T R - QB 


(3. 413CALO)! + (0. 09459 + 3. 413CDLO)! + 13. 3(1 - E cg ) + QA 


(147) 


For the transient condition, by substituting equations (146) and (112) into equation (70), 
the governing equation takes the form of equation (75a) with the constants A, B, and 
C given by equations (122), (123), and (124); however, in equations (121), (122), and 
(123) the voltage curve fit constants CAHI, CBHI, CDHI, and CEHI are replaced by 
CALO, CBLO, CDLO, and CELO, respectively. The value of F for case IV is the 
same as in equation (146). The solution to the differential equation is given by equa- 
tion (82), with the value of T the same as in equation (147). 

s s 


Case V 

The conditions for case V are as follows: 

1. 0<is5 amperes 

9 t < T < T 

i BPLO _ ~ BPHI 

3. Partial regeneration ( <p £ 0) 

4, In-line heater off 

For case V, the only three parameters which are different from case II are £ 3 * 

and 6 ^. These are derived using the same techniques as were used in case n, except 

that equation (146) is used for the voltage instead of equation (119). The equations for 
S 2 , £ 3 , and 6 4 are 


i 2 = (3. 413CBLO)I 2 + (3. 413CDLO + 0. 01043)1 + 6 2 + QA (148) 

l 3 = (3. 413CBLO)! 2 + (3. 413CELO - 133. 3)1 + 6 g + QB (149) 

6 4 = -(3. 413CALO)! 2 - (3. 413CDLO + 0. 01043)1 - QA (150) 
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and the value of r is the same as that used in equation (146). The steady -state and 
transient temperature equations are the same as for case II, except that the new values 
of £ 2 , and r are used for case V. 

The foregoing equations were used to formulate a computer program which ac - 
cepts a given current profile and predicts the thermal response of the fuel cell to the 
demand load imposed upon it. The computer program is described in appendix B. A 
computer program was also written for the simplified analysis presented in the section 
on simplified transient analysis. 


COMPARISON OF ANALYTICAL RESULTS 


Before comparing the results of the general analysis to the test data obtained in 
the thermal vacuum chamber, the results of the simplified analysis and the general 
analysis will be compared. Figure 21 illustrates this comparison for a sample current 
profile. It is evident from the figure that the incorporation of the primary regenerator 
and the bypass valve into the analysis, which is the major difference between the sim- 
plified and general analyses, has two main effects. First, by using the equations of the 
general analysis, the thermal response time of the fuel cell is much smaller in the re- 
gion of regeneration than the response time obtained from the simplified analysis. Sec- 
ond, by using the equations of the general analysis, the operating -temperature range of 
the fuel cell is narrowed considerably, particularly at high -power loads. 

The steady -state voltage -current 
characteristic of the fuel -cell module 
(obtained by using the equations of the 
general analysis) is compared in figure 22 
to the thermal -vacuum characteristics 
(obtained from test data taken from in- 
house performance evaluations of Apollo 
fuel -cell power systems) of two qualified 
Apollo fuel -cell modules. The steady - 
state voltage -current characteristic ob- 
tained with the simplified analysis is also 
included for comparison. In this figure, 
the large error incurred by use of the 
simplified analysis is shown clearly. As 
is evident in the figure, the general anal- 
ysis gives a more realistic approximation 
to actual fuel-cell performance. 


Figure 21. - Comparison of results of 
the simplified analysis and the 
general analysis. 


The empirical data used in the equa- 
tions are. a statistical average of data ob- 
tained from many actual power plants and 
components. Because each power plant 
is a separate entity, the performance of 


NASA-S-67-5551 
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each will be slightly different. Therefore, 
the mathematical model can only be ex- 
pected to predict average power plant per- 
formance. Examples of performance 
differences between actual power plants 
are indicated by variations in the steady - 
state voltage -current characteristic 
(fig. 22). 

The discontinuity in the characteris- 
tic curves (fig. 22) of the general analysis 
and the qualified power plants is that point 
at which the primary bypass valve is fully 
closed and the power plant is in a condition 
of full regeneration. If the current is fur- 
ther decreased from this point, no more 
regeneration is obtained, and the stack 
temperature settles to a lower steady - 
state value, thereby lowering the steady- 
state voltage. 

The constant -temperature lines 
shown in the voltage -current plot in fig- 
ure A-l of appendix A for a constant 
electrolyte (P^O) concentration of 27 per- 
cent are not included in figure 22 because 
at each different point on the steady -state 
characteristic line, the electrolyte concentration has a different value. If lines of 
constant temperature and constant concentration were both shown, a three-dimensional 
plot would be needed. 

In figures 23(a) to 23(d), the transient performance results of the general analy- 
sis are compared to actual test data obtained for various step load changes from in- 
house performance evaluations of a qualified Apollo fuel-cell module. Note in making 
the comparison, that the steady- state voltage values of the general analysis will be 
slightly higher in the regeneration region than those of the qualified power plant to 
which they are compared in figure 22. At low currents, however, the two steady-state 
lines approach each other. 

In the figures, the glycol-coolant inlet temperature to the module and the glycol 
flow rate through the condenser are given for the test module. These two quantities 
determine the primary-side condenser exit temperature for any given fuel-cell 

load. The maximum and minimum values of T c / g observed during the various test 

sequences are also shown on the figures for comparison with the constant condenser 
exit temperature of 160 ° F assumed in the general analysis. The variations in con- 
denser exit temperature are partially caused by the fact that as fuel-cell load increases 
more water is produced, with a resultant increase in the partial pressure of the water 
vapor in the primary loop and a consequent increase in condenser exit temperature. 
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Figure 22. - Comparison of steady-state 
voltage current characteristics. 
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Figure 23(a) indicates a step from 30. 0 to 34. 5 amperes with condenser exit tem- 
perature varying between 159. 9° and 166. 3° F on the test module. The steady-state 
voltage of the test module is approximately 0. 1 volt lower than that of the analytical 
model, while the steady -state temperature of the test module is higher than that of the 
analytical model by about 5° F. As illustrated in figure 23(b), the stack temperature 
was allowed to settle for 2 hours at 45 amperes; then the power plant was put on open 
circuit. The analytical model and the test module compared favorably in temperature 
just before the load was increased to 45 amperes. Six hours later, the load was de- 
creased to 37 amperes. This also is illustrated in figure 23(b). Note that figure 23(c) 
also illustrates the temperature and voltage response comparison for a step load de- 
crease in current. 

The first part of figure 23(d) illustrates the operation of the low-power in-line 
heater. The temperature dead band of the analytical model was set between 385° and 
387° F and thus compared favorably with the test module. Note, however, that the 
voltage drop in the analytical model when the heater is energized exceeds that of the 
test module. This is caused partly by the use of a smaller heating -element resistance 
in the analytical model, which causes the model to draw a larger heater current at 
some nominal voltage and results in a larger voltage drop. The larger voltage drop in 
the analytical model is also caused by the selection of a nominal module voltage of 
31 volts (assumed constant) when the in-line heater is energized. The in-line heater of 
the analytical model was energized more frequently than that of the test module because 
the fuel-cell environment temperature of the test module was 150° F; whereas, an en- 
vironment temperature of 77° F was used in the general analysis. This increased en- 
vironment temperature of the test module caused a lower stack heat loss than the 
environment temperature used in the general analysis. 

Figure 23(d) illustrates a very interesting feature of the general analysis. When 
the load is increased from 18 to 25 amperes at t = 4 hours, the stack temperature for 
the analytical model has just dropped to 385° F, which causes the heaters to be ener- 
gized with simultaneous application of the step load. This results in a total instantane- 
ous current demand of 12. 15 amperes (since for the analytical model the heater current 
is assumed constant at 5. 15 amperes), thus causing a voltage drop to 27. 1 volts. In 
the test module, however, the in-line heater was still off when the step load was applied, 
and thus the voltage only dropped to 28. 6 volts. Several minutes after t = 4 hours, the 
stack temperature of the analytical model reached 387° F, which caused the heater to 
be deenergized and the voltage to rise instantaneously to 28. 55 volts, from which point 
the normal transient voltage buildup continued. At t = 8 hours, the power plant was put 
on open circuit for 31 minutes. After this time, the 25-ampere load was resumed. 


CONCLUDING REMARKS 


From the steady-state voltage-current-characteristic comparison, it is evident 
that the general analysis yields a characteristic curve which compares favorably with 
the curve for the qualified Apollo fuel -cell module used for transient comparison. The 
comparison of the transient response of the general analysis to the response of the qual- 
ified fuel -cell module yields the conclusion that the equations derived in the general anal- 
ysis and the representative performance characteristics which served as inputs to the 
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analysis are adequate to describe the performance of an actual fuel -cell system through- 
out its normal operating range for both the steady -state and the transient modes of op- 
eration. 

The following is a summary of the principal limitations of the analysis. 

1. The predicted values of the analysis are only valid for a fuel-cell power plant 
operating in a vacuum environment. 

2. The analysis is restricted to a normal range of glycol -coolant inlet temper- 
atures such that the condenser exit temperature on the primary side remains within 
approximately ±10° from 160° F. 

3. The environment temperature of the fuel cell is restricted to temperatures in 
the region of 77° F. It is thought that a variation of ±20° F from 77° F would be toler- 
able. 


4. The stack heat loss is assumed to have only a linear variation with temper- 
ature. Little is known about the stack heat loss, since the stack and the primary 
regenerator are inside the pressure vessel where it is very difficult to place instru- 
mentation. However, it is thought, and test -data comparisons seem to agree, that the 
linear assumption is adequate for a temperature of 77° F. 

5. The cold-side effectiveness of the primary regenerator was assumed to be 
constant at 84 percent. (Cold-side effectiveness is known to be a function of temper- 
ature, percent bypass, and specific humidity. ) This limitation and limitation 4 do not 
appear to cause serious errors in the computations, however. 

6. By assuming a purely resistive load, Ohm's law dictates that the current 
through the load is a function of the voltage impressed upon it which, in the actual 
spacecraft, is fuel-cell output voltage. However, the fuel-cell output voltage changes 
with fuel-cell temperature, fuel-cell concentration, and other variables, as observed 
in the analysis, which causes fluctuations in fuel-cell output current. Although these 
fluctuations are usually small, this effect should be taken into account. This was not 
done in this analysis. 

In conclusion, it is believed that this analysis presents an adequate mathematical 
model for predicting fuel-cell performance, subject to the limitations imposed upon the 
study. If performance predictions are desired for power plants other than the fuel-cell 
used for transient comparison, the performance characteristics of the specific power 
plant under consideration must serve as an input to the computer program. 

Areas of further study in order of priority are as follows: 

1. Taking the effects of variable current into consideration would mean an 
iterative -type solution in which, at any given time, a fuel-cell output voltage is as- 
sumed for a given resistance input. From this assumption a current is calculated. By 
using this current as an input to the general analysis, fuel-cell temperature and voltage 
are calculated for the given time. The calculated voltage is compared to the assumed 
voltage. If the two are not equal, some new voltage is assumed, and the process is 



repeated until the calculated and assumed voltages are within a reasonable tolerance. 
The computation then proceeds to the next time step, and the whole iteration process 
is repeated. 

2. Taking into account the variations of fuel -cell performance with changes in 
environment temperature would include trying to define the stack heat loss more ac- 
curately, since stack heat loss is a function of environment temperature and stack 
temperature. 

3. The effect of voltage degradation with time caused by the buildup of impu- 
rities within the cells is another area of further study. 

4. Taking into account the variations in fuel -cell performance caused by changes 
in condenser exit temperature on the primary side would entail an analysis of con- 
denser and secondary -loop performance. 

The analytical techniques required to formulate a mathematical model for a 
complex system such as the fuel -cell system may at first appear to be unique and, 
therefore, lacking in general applicability; however, there is a common body of knowl- 
edge in modeling and simulation which is broadly applicable and of great help to the 
systems engineer. This report is only one application of these general techniques, but 
the approach taken and the methods used are sufficiently general to allow application to 
other types of systems. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 12, 1968 
941-13-05-01-72 
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APPENDIX A 


DEVELOPMENT OF EMPIRICAL RELATION FOR FUEL-CELL 
POWER PLANT OUTPUT VOLTAGE AS A FUNCTION OF CURRENT, 
TEMPERATURE, AND ELECTROLYTE CONCENTRATION 


To obtain the necessary relation for the module output voltage, a statistical com- 
pilation of test data (triangular points) taken from many power plants is plotted 
(fig. A-l). Test data are not readily available for the open-circuit condition at tem- 
peratures above 440° F or for the current region between open circuit and approxi- 
mately 15 amperes. Therefore, values in this region can only be estimated. The 
comparison in figure A-l is made for a constant electrolyte concentration of 27 percent 
water by weight. 
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Figure A-l. - Comparison of polarization curves of analytical 
model with those of statistical test data. 
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Since the fuel cell normally operates at temperatures between 380° and 440° F, 
at currents between 20 and 55 amperes, and at voltages between 27 and 31 volts, only 
the data points in the shaded area in figure A- 1 were considered for the current region 
from 5 to 55 amperes. From these points, it is evident that at any specified tempera- 
ture, the module output voltage is essentially a linear function of current. The voltage 
equation would thus have the form 


E = nijl + b^ 


(Al) 


where m^ is the slope, and b^ is the Y-intercept. 

By assuming linearity, the slope and the Y-intercept are calculated for each 
constant-temperature line between 380° and 440° F. Figure A-2 is a plot of m^ 

versus T and b^ versus T. Calculations indicate that the use of a first-order 

approximation in the most-frequented operating-temperature range (400° to 440° F) 
produces only a very small error in the voltage calculation for the region between 380° 
and 400° F. This approximation yields relations for the slope and Y-intercept as a 


function of temperature as follows. 




Temperature/ °F 

Figure A-2. - Variation of slope and 
Y-intercept of constant-temperature 
voltage lines with temperature. 


m 1= 7.96 X 10" 4 T - 0.481 (A2) 


b x = 0. 026T + 23. 0 (A3) 


By substituting equations (A2) and (A3) into 
equation (Al), the voltage variation E y 

with current and temperature at a constant 
electrolyte concentration of 27 percent 
water is 


E y = (7. 96 X 10" 4 T - 0. 48 l)l 

+ (0. 026T + 23. 0) (A4) 


where 5 amperes < I < 55 amperes, and 
380° F < T 440° F. 

It is evident in figure A- 1 that agree- 
ment between the test data and equation (A4) 
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is very good within the region considered. There is also good agreement considerably 
outside of this region. 

To account for changes in electrolyte concentration, a correction factor K c is 
introduced. Thus, the general voltage equation can be written as 

E = E v + K £ (percent HgO - 27. 0) (A 5) 


Extensive tests conducted on single cells, on six-cell stacks, and on complete modules 
have shown that K c is essentially constant over the normal electrolyte-concentration 

operating range of the fuel cell, which is less than ± 2 percent HgO about the nominal 

27 percent H 0 0. The value of K from these tests is -0. 2473 V/A percent H-O above 

27 percent. Thus, the general equation for the variation of module output voltage with 
current, temperature, and electrolyte concentration is 


E = (7. 96 X 10" 4 T - 0. 48 l)l + (0. 026T + 23. 0) 
- 0. 2473 (percent H^O - 27. 0) 


(A6) 


where 5 amperes < I < 55 amperes, 380° F s T ss 440° F, and 22 s percent s 32 . 

It is evident from equation (A6) that an increase in fuel-cell operating (stack) tempera- 
ture will raise the voltage output of the module correspondingly. For the temperature 

-4 

range considered, the term 7. 96 x 10 T - 0.481 is always negative; therefore, a 
current increase causes a decrease in module output voltage. 

The quantity percent HgO is calculated using the assumption of equilibrium for 

the partial pressure of the water vapor above the KOH within the stack. Under this 
assumption, the electrolyte concentration is a function of (1) the partial pressure of the 
water above the KOH and (2) the temperature. Figure A- 3 is a plot of a statistical av- 
erage for percent KOH as a function of these two variables taken from actual laboratory 
data (triangles). For a nominal fuel-cell condenser exit temperature of 160° F on the 
primary side, the partial pressure of the water above the KOH will always be between 
4. 6 and 7. 0 psia. A first-order approximation biased to the most-frequented region is 
again used. Figure A-4 is a plot of the variations of slope and Y-intercept. The slope 
is approximately constant at -1.0 percent change in KOH per unit change in partial 
pressure. By using a straight line fit for the Y-intercept in figure A-4, the relationship 
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Partial pressure of water vapor above KOH, psia 


Figure A-3. - Variation of electrolyte concentration 
with temperature and partial pressure. 
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Figure A-4. - Variation of Y-intercept of constant-temperature 
concentration lines with temperature. 


for electrolyte concentration as a function of partial pressure and temperature is 
percent KOH = 0. 0997T + 38. 22 - p , from which 


percent HgO = 61. 78 - 0. 0997T + p^ 


(A7) 


where 350° F ^ T ^ 500° F, and 4. 6 psia < p^ < 7. 0 psia. Equation (A7) is used in 
conjunction with equation (A6) in the analysis. 

Since very little data could be obtained for this analysis between open circuit and 
20 amperes, the linear fit in equation (A4) is used down to 5 amperes. Then, using 
test-data values for open-circuit voltage, the same procedure used in the derivation 
of equation (A4) is used for the low-current region (0 to 5 amperes), with the result 
that 


E = (2. 775 X 10" 3 T - 1. 903)l + (0. 0157T + 30. 24) 
- 0. 2473 (percent H^O - 27. 0) 


(A8) 


where the same correction factor used before for concentration changes is applied. 
Percent H 2 0 is calculated from equation (A7), the limits of which are sufficient to in- 
clude operation down to open-circuit conditions. The solid lines in figure A-l are a 
plot of E from equation (A8) for various temperatures. 

Equations (A6) and (A8) represent the empirical relations used in the analysis for 
the high-current (5 to 55 amperes) and low-current (0 to 5 amperes) regions, respec- 
tively. These two equations should be used with equation (A7) for percent K^O. The 

quantity p w is calculated from equation (56). 
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APPENDIX B 


THE COMPUTER PROGRAM FOR THE GENERAL ANALYSIS 


A digital computer program was written for the simplified analysis presented in 
the section on simplified transient analysis and was later expanded for the general 
analysis presented in the section on general transient analysis. The program for the 
general analysis is presented in this appendix. 

The digital program accepts a current-time profile and predicts for this profile 
the transient temperature and voltage response of the fuel-cell module. The program 
includes the five current and temperature regions of operation for which transient 
equations were developed in the section on general transient analysis. 


SYMBOLS 


AGAMA 


collection of terms used in the equations for simplifica- 
tion 

collection of terms used in the equation for simplification. 
AGAMA is a function of GAMA, and is denoted by r 
in the analysis 

collection of terms used in the equations for simplifica- 
tion 


BP 

BP1, BP2 


percent bypass around primary regenerator 

curve fit constants for primary-bypass-valve character- 
istic 


collection of terms used in the equations for simplifica- 
tion 


Cl 


collection of terms used in the equations for simplifica- 
tion 


CAHI, CBHI, CDHI, CEHI 


curve fit constants for fuel-cell performance character- 
istics in the high-current region, 5 A < I 55 A 


CALO, CBLO, CDLO, CELO curve fit constants for fuel-cell performance character- 
istics in the low- current region, 0<I<5A 


COLD 


Wc. 


for hydrogen and water vapor at cold side to 


primary regenerator, Btu/° F 
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DELHT 


DISCRM 

DT 

DTI 

E 

EPR 

ESS 

G 

GAMA 


HTIME 

I 

II 
12 

IHTR 

LAMBDA 

PCH20 

PSS 

PW 

Q 

QA, QB 
RT 
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length of time the in-line heater stays on, hr 
discriminant of the polynomial in W^, and 


time increment for transient temperature and voltage 
calculations, hr 

dummy variable (time) increment used in calculating the 
time at which the bypass valve closes, hr 

module voltage, V 

primary regenerator cold- side effectiveness (fractional) 

steady- state module voltage, V 

dummy variable, equal to y/q or v / c q, depending on 
whether q is positive or negative 


specific humidity of hydrogen and water-vapor mixture 


at stack outlet, 


lb H 0 0 
m 2 

lb dry H 
m J 


2 


time index used to calculate the length of time the in-line 
heater stays on, hr 

total fuel-cell current after a given step load change 

P'V'htr)- A 

current before step load change, A 

current after step load change, A 

current drawn by in-line (low-power) heater, A 


time constant for fuel-cell transient response, hr 
percent H^O in electrolyte solution 


steady- state power level, W 

partial pressure of water vapor at stack outlet, psia 

negative discriminant of polynomial in A, B, and C 

curve fit constants for stack heat loss 

variable time used to calculate the transient temperature 
and voltage values between TIME1 and TIME2, hr 



RT1 

S 

T 

T1 

TAN 

TB 

TBPHE 

TBPLO 

THETA, THETA1, THETA2 

THOF 

THON 

TIME 

TIME1 

TIME 2 

TO 

TR 

TSS 

TSTART 
Wl, W2, W3 


dummy variable (time) used in calculating the time at 
which the bypass valve closes, hr 

time index used in calculating the time at which the pri- 
mary bypass valve closes, hr 

fuel-cell (stack) temperature at any time, °F 

initial fuel-cell temperature at each step load change, °F 

collection of terms used in the equations for simplifica- 
tion 

dummy variable (temperature) 
temperature of bypass, high 
temperature of bypass, low 

collection of terms used in the equations for simplifica- 
tion 

temperature at which the in-line heater is deenergized, 

°F 

temperature at which the in-line heater is energized, °F 

dummy variable used in calculating the time at which the 
primary bypass valve opens, hr 

dummy variable used to initialize time at each step load 
change, hr 

time at which the load changes from 12 (or I, if in-line 
heater is on) to some new value of current (12 equals 
fuel-cell current from TIMEl to TIME2), hr 

initial time of current 12 (or I, if in-line heater is on) 
(TO = TIMEl), hr 

temperature of hydrogen and water-vapor mixture at 
cold- side inlet of primary regenerator, °F 

steady-state fuel-cell (stack) temperature, °F 

initial temperature of fuel-cell stack at time t = 0, °F 

collection of terms used in the equations for simplifica- 
tion; denoted by (jp an d £3 in the analysis 
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XI, X2, X3, X4 
Yl, Y2, Y3 
Z 

Zl, Z2, Z3, Z4 


collection of terms used in the equations for simplifica- 
tion 

collection of terms used in the equations for simplifica- 
tion; denoted by /3j, /Sg, and /3g in the analysis 

term used in the equations for simplification; denoted by 
a in the analysis 

collection of terms used in the equations for simplifica- 
tion; denoted by 6j, Sg, and 5^ in the analysis 
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PROGRAM MAIN 

REAL II. 12. IHTR. I .LAMBDA 

WRITE16.1) 

1 FORMAT ( 55H1TRANSIENT THERMODYNAMIC ANALYSIS OF A FUEL CELL SY! 
17H SIMON //) 

50 READ (5.51) CAHI »CBHI »CDHl .CEHI . C ALO » CBLO ► CDLO . CELO 

51 FORMAT ( 8F10 . 0 ) 

IF(CAHI.EQ.O.O) GO TO 2000 
READ ( 5 . 51 ) BP1 * BP2 . TBPLO . TBPHI 
READ (5.9) EPR 

READ ( 5 . 51 ) QA » QB 
READ (5.51) THON » THOF 
70 READ (5.9) DT 
READ (5»9)TSTART 
9 FORMAT(FIO.O) 

TIME1=0.0 
T1=TSTART 
T-T 1 

5 READ (5.10) TIME2.I2 
10 FORMAT(2F10.0) 

IF(TlME2)bO» 70.80 
80 WRITE ( 6 » 999 ) 

T0=TIME1 
R f=T0 

2 IF ( T . LT • 350 • ) GO TO 24 

TK= ( 0 . 881 * 12 + 2261 . ) / ( 8. 84 IE-O 3 * 12+13. 3) 

GAMA=(2.70+22.976E-03*I2)/3.5l 
Pw=60. 0/(1. 0+8. 94/GAMA) 

AGAMA=-8. 601-14. 838/ <1.0+8. 94/GAMA) 

IFCT.LE. TBPLO) GO TO 3 
IHTR=0.0 

17 IF(T. GT. TBPHI ) GO TO 4 
C0LD=8. 84 IE-03+12+1 3. 3 
Z=10.431E-03*I2+13.3 

yi=-epr*z*bpi*.oi/cold 

Y2=(EPR*Z/COLD)*(1.0-BP2*.01+BP1*TR*.01) 
Y3=(EPR*Z*TR/COLD>*(BP2*. 01-1.0) 

Zi=-13.3*Y1 

Z2=13.3*(1.0-Y2) 

Z3=-13.3* (TR+Y3) 

W1=Z1 

IF (12. LT. 5.0) GO TO 34 

W2=Z2+Q.08416*I2+3.413*I2*CDHI+3.413*I2**2*CAHI+0.01043*I2+QA- 
W3=Z3 +3.413*I2*AGAMA+3.413*I2*CEHI+3.413*I2**2*CBHI-13. 

12+QB 

28 DISCRM=W2**2-4.0*W1*W3 

IF (DlSCRM.LT .0.0) GO TO 20 
TSS=(-W2+SQRT(DISCRM) ) / ( 2 • 0*W1 ) 

PCH20=61.78-G.0997*TSS+PW — ~ - - — - - 

IF( I2.LT.5.0) GO TO 38 

ESS= ( CAHI*TSS+CBHI ) *12+ ( CDHI*TSS+CEHI ) -0.2473* (PCH20-27.0) 

GO TO 39 

38 ESS=(CALO*TSS+CBLO)*I2+(CDLO*TSS+CELO)-0.2473*(PCH2O-27.0> 

39 IF( I2.GE.5.0) GO TO 35 

Z4=-3. 413*CAL0*I2**2+(-3.413*CDL0-0. 01043-0. 084161*12 - 

GO TO 29 
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35 Z4=-3.413*C AH I *12**2+ (-3. 413*CDHI-0 . 0l0<+3-0 . 08416) *12 
29 C1=Z1*TSS**2+(Z2-Z4+QA)*TSS 
A=C1 

B=Z4-Z2-QA 

C=-Z1 

Q=4.0*A*C-B**2 
If (Q)6»7»8 
6 G=SQRT(-Q) 

X1=<2.0*C*T1+B)/G 
14 X2=-TANH( lG*(RT-T0) J/64.0) 

T=-B/(2.0*C)+(G/(2.0*C) )*( ( X1+X2 ) / ( 1 . 0+Xl*X2 ) ) 

PCH20=61.78-0.0997*T+PW 

IF ( 12. LT .5.0) GO TO 4q 

£-(CAHI*T+CBHI ) *12+ (CDHI*T+CEHI ) -0 . 2473* ( PCH20-27 . 0 ) 

GO TO 11 

40 E=(CAI_0*T+CBLO)*I2+(CDLO*T + CELO)-0.2473*(PCH20-27.0) 

11 IF (T.GT.TBPHI )GO TO 4 
IF(T.L£.TBPLO)GO TO 12 
BP=BP1*T+BP2 

WHITE (6. 1000) FT. l2»E*T»bP» PCH20 » IhTR 
IF (RT.GE.TIME2)G0 TO 60 
R|=RT+DT 
GO TO 14 

34 W2=Z2+0 . 08416*12 + 3. 413*I2*CD|_0+3. 413* 1 2**2*C ALO+O . 0 1 043*12+0 A 
W3=Z3 +3.413*I2*AGAMA+3.413*I2*CEL0+3.413*I2**2*CBL0-13. 

12 + QB 
GO TO 28 
4 WRITE (6» 2001 ) 

2001 FORMAT ( /// 15X28H STACK TEMP. EXCEEDS 450 F ) 

GO TO 2000 

12 IF(Q.GT.O.) GO TO 13 
RT=RT-3T 

S=RT 

DT1=DT/10.0 

RT1=RT+DT1 

21 X2=-TANH( (G*(RT1-T0) )/64.0) 

T--B/ (2.0*C)+(G/(2.0*C) )*( (X1+X2) / C1.0+Xl*X2) ) 

RT=RT1 

xxxx=tbplo-.ooi 

IF (T.GT.XXXX) GO TO 22 

Rr=RT-DTl 

DT1=DT1/10.0 

RT1=RT+DT1 

GO TO 21 

22 YYYY=TBPLO+.001 
IF(T.LT.YYYY) GO TO 15 
RT1=RT1+DT1 

GO TO 21 

15 PCH20=61.78-0.0997*T+PW 
IF (I2.LT.5.0) GO TO 41 

E=(CAHI*T+CBHI ) *I2+(CDHI*T+CEhI ) -0.2473*(PCH20-27.0) 

GO TO 42 

41 E=(CALO*T+CBLO) *12+ (CDLO*T+CELO 5-0.2473* (PCH20-27.0) 

42 Ti=T 
BP=0.0 
IHTR=0.0 
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WRITE (6 » 1000 )RT1 , I2r£»T»BP'PCH20*IHTR 

IF (RT1.GE.TIME2) GO TO 60 

RT=S+DT 

Tl=TBPLO 

T0=RT1 

IF (T.LE.THON)GO TO 16 
IHTR=0.0 

IF( I2.GE.5.0) GO TO 25 

TR=( 0.681* I 2+2261. ) / ( 8 . 841E-03*I2+13 . 3 ) 

GAMA= (2.70+22. 976E-03*I2 ) /3 • 51 
Pw=60./(1. +8. 94/GAMA) 

AGAMA=-8.601-14.638/ ( 1 . 0+8 . 94/GAMA ) 

TbS = ( -3 . 413*CBL0*I 2**2+ (133.3-3 . 4l3*CEL0 — 3 . 413* AG AM A )*I2+13.0 : 
1-EPR) *TR-QB>/( 3.41 3*CAL0* 12**2+ < • 09459+3 . 4 13*CDL0 ) *12+13 . 0* ( 1 
2R ) +QA ) 

PCH20=61 . 78-0 . 0997*TSS+PW 

EbS=(CALO*TSS+CBLO)*I2+(CDLO*TSS+CELO)-0.2473*(PCH2O-27.0) 

A=32.0 

B = 13.3*( 1.0-LPR) + ( 3.41 3*CALO) *12**2+ (3.41 3*CDL0+3.'+l 3*0.0247 ) : 
1A 

LAMBDA=B/A 

27 T=TSS-(TSS-T1)*EXP(-LAMBDA*(RT-T0) ) 

PCH20=bl .78-0.0997*T+PW 

E-(CAL0*T+C6L0)*I2+(CDL0*T +CELO )~0.2473*( PCH20-27 . 0 ) 

BP=0.0 

IF(T.GT.TBPLO) GO TO 32 

IF (T.LE.THON) GO TO 33 

I F ( T . LT . 350 . ) GO TO 24 

WRITE (6 r 1000 )RT. I2»Ef T»BPrPCH20 » IhTR 

IF (RT.GE.TIME2) GO TO 60 

RT=RT+DT 

GO TO 27 

2b TR=( 0.681 *12+2261. ) / ( 8 . 841E-03* I 2+13 . 3 ) 
GAMA=(2.70+22.976E-03*I2) /3.51 
PW=60.0/ ( 1 .0+8.94/GAMA) 

AGAMA-~8.601~14.838/ ( 1 . 0+8 . 94/GAMA > 

TbS=(-3.413*CBHI*l2**2+(133.3-3.+13*CEHI-3.4l3*AGAMA)*I2+13.0^ 
1-h.PR) *TR“QB)/(3.413*CAHI* 12**2+ ( . 09459+3 . 413*CDH I ) *12+13 . 0* ( 1 
2R)+QA) 

PCH20=61 .78-0 . 0997*TSS+PW 

EbS=(CAHI*TSS+CBHI ) *12+ ( CDHI*TSS+CEHI ) -0 . 2473* ( PCH20-27 . 0 ) 
A=32.0 

B=13.3*(1.0-EPR)+(3.413*CAHI)*I2**2+(3.413*CDHI+0.09459)*I2+O, 

LAMBDA=B/A 

18 T=T^S-(TSS-T1)*EXP(-LAMbDA*(RT-T0) ) 

PCH20=61 . 78-0 . 0997*T+PW 

E=(CAHI*T+CBHI ) *12+ (CDHI*T+CEHl) -0.2473* (PCH20-27.0) 

BP=0.0 

IF (T.GT.TBPLO) GO TO 37 
IF(T.LE.THON) GO TO 36 
WRITE ( 1000) RTr I2.E»TfBP»PCH20» IHTR 
IF ( RT • GE • TIME2 ) GO TO 60 
RT=RT+DT 
GO TO 18 
37 TB=TBPLO 
IhTR=0.0 
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60 TO 19 
36 Tb=THON 
IHTR=5.1b 
GO TO 19 

33 TIME=T0“ ( 1 ./LAMBDA ) *ALOG ( ( T5S-TH0U ) / ( TSS-T1 ) ) 

T=TH0N 

PCH20=61 . 78-0 . 0997+T+PW 

IHTR=5.1b 

I=l2+IhTR 

E=(CALO*T+CBLO)*l2+(CDLO*T+CELO)-U.2473*(PCH2O-27.0) 

BP=0.0 

WRITE (6* 1000) TIME* 12* E*T* BP* PCH20, IHTR 
E=(CALO*T+CBLO)*I + ( CDLO*T+CELO ) -0 .2473* ( PCH20-27 . 0 ) 

WRITE (6* 1000) TIME* I * E * T * BP * PCH20 * IHTR 

T1=T 

T0=TIME 

HTIME=TIME 

Gu TO 26 

32 T1ME=T0- 11. /LAMBDA ) *ALOG ( l TSS-TBP lO ) / ( TSS-T1 ) ) 

T=TBPLO 

PCH20=61.78-0.0997*T+PW 

E=(CAL0*T+C8l0> *I2+< CDLO*T+CELO> -0.2473* (PCH20-27.0) 

BP=0.0 

WRITE <6* 1000) TIME* 12 * E * T . BP * PCH20 , IHTR 

T1=T 

T0=T I ME 

IF ( TIME »GE. TIME2 ) GO TO 60 
GO TO 17 
16 IHTR=b.lb 
I=l2+b.lb 

TR= (0.681* 1+2261. ) /(8.841E-03*I+13.3) 

GAMA=(2.70+22.976E-03*I)/3.51 

PW=60./ Cl. +8. 94/GAMA) 

AGAMA=-8. 60 1-14. 838/ (1.0+8. 94/GAMA) 

TSS=< 13.3* (1.0-EPR)*TR+544. 8-QB+133. 3*1-3. 413*1* (CBHI*I+CEHI+, 
1) )/(13.3* t 1.0-EPR) +0. 01043*1 +3. 413* I* ( C AHI*I+CDHI+. 02466 )+GA) 
PCH20=6l • 78-0 . 0997*TSS+PW 

ESS=(CAHI*TSS+CBHI ) *1+ ( CDHI*TSS+CEHI ) -0 . 2473* ( PCH20-27 . 0 ) 
A=32.0 

B=13.3*(1.0-EPR)+3.413*CAHI*I**2+(3.413*CDHI+0.09459)*I+QA 

LAMBDA-B/ A 

DELHT=RT-TIME 

23 T=TSS- ( TSS-T 1 ) *EXP ( -LAMBDA* ( RT-TO ) ) 

PCH20=61.78-0.0997*T+PW 

E= ( CAHI *T+CBHI ) * 1+ (CDHI*T+CEHI ) -0.2473* ( PCH20-27 . 0 > 

BP=0.U 

IF(T.GT.THOF) GO TO 30 

IF ( T . LT • 3b0 . 0 ) GO TO 24 

WRITE (6* 10 00) RT* 12 * E * T * BP * PCH20 * IHTR 

IF (RT.GE.TIME2) GO TO 31 

RT=RT+DT 

DELHT=DELHT+DT 

GO TO 23 

31 WRITE ( 6 * 1004 ) DELHT 
GO TO 60 

30 TIME=T0“ ( 1 • 0 /LAMBDA ) *ALOG ( ( TSS-THOF )/( TSS-T 1 ) ) 


78 



FTN 5.3 


05/05/67 


D£LHT=DELHT-(RT-TIME) - - — 

T=THOF 

PCH20=61 . 78-0 . 0997*T+PW 

IHTR=5.15 

I=I2+IHTR 

E=(CAHI*T+CBHI)*I +(CDHI*T+CEHl)-0.2473*(PCH2O-27.0) 

BP=0.0 - - 

WRITE (6» 1000) TIME . I »E»T .BP.PCH20. IHTR 
E=(CAHI*T+CBHI)*I2+(CDHI*T+CEHI)-0.2473*(PCH20“27.0) 

WRITE (6. 1000 ) TIME. I2.E»T»BP*PCH20.IHTR 

WRITE(6»1004) OELHT 

T0=TIME 

26 T1=T - -- 

IF(TIME.GE.TIME2) GO TO 60 
GO TO 3 

19 TIME=T0- ( 1 . /LAMBDA ) *ALOG ( (TSS-TB ) / ( TSS-T1 ) ) 

T=TB 

PCH20=61 • 78-0 . 0997*T+PW 

E=(CAHI*T+CBHI ) *I2+(CDHI*T+-CEHI >-0.2473* (PCH20-27.0) - ' 

BP=0.0 

WRITE (6. 1000) TIME. 12 »E»T.BP»PCH20. IHTR 

T1=T 

T0=TIME 

IF(TIME.GE.TIME2) GO TO 60 

IF ( T . LE .400 • ) GO TO 16 - - - — -- 

GO TO 17 

20 WRITE ( 6 . 2003 ) 

2003 FORMAT (///15X31H DISCRIMINANT IS LESS THAN ZERO) 

GO TO 2000 

7 WRITE ( 6 » 2004 ) 

2004 FORMAT ( ///16X25H DISCRIMINANT EQUALS ZERO) - - - 

GO TO 2000 

13 X3=ATANl ( <2.0*C*420.0)+B)/G> 

X4=ATAN( ( <2.0*C*T1)+B)/G) 

TIME=T0+(64.0/G) *(X3-X4) 

T=TBPL0 

RT1=TIME - — - 

GO TO 15 

8 THETA1=ATAN( ( 2 . 0*C*T1+B ) /SORT ( Q ) ) 

THETA2=(SQRT(Q)*(RT-T0) J/64.0 
THETA=THETA1+THETA2 

TAN=SIN ( THETA) /COS ( THETA ) 

T=(-B)/(2.0*C)+(SQRT(Q)/(2.0*C) >*TAN - — 

PCH20=61 . 78-0 . 0997*T+PW 
IF ( I2.LT.5.0) GO TO 43 

E= (CAHI*T+CBHI ) *12+ (CDHI*T+CEHI )-0 .2473* (PCH20-27. 0) 

GO TO 44 

43 E=(CALO*T+CBLO) *12+ (CDLO*T+CELO) -0.2473* (PCH20-27.0) 

49 G=SQRT(Q) 

GO TO 11 

24 WRITE (6.2005) 

2005 FORMAT (///15X27H STACK TEMP. IS BELOW 350 F) 

GO TO 2000 

60 PSS=ESS*I2 
T1ME1=TIME2 
11 = 12 
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T1=T 

WRlTE(6»1001)TSS 
WRITE 1002) ESS 
WRITE(6.1003)PSS 
GO TO 5 

999 FORMAT (7X4HTIME7X1HI11X1HE11X1HT8X2HBP5X4HC0NC5X4HIHTR/7X3HHR: 
1AMPS7X5HV0LTS7X4HDEGF6X5HPRCNT2X7HPRCTH203X4HAMPS) 

5000 FORMAT (//4F15.7) 

1000 FORMAT (F12.4»2FH.4»F12.4»F8.2*F9.2»f8»2) 

1001 FORMAT (//14X20H TEMP STEADY STATE = F13.8) 

1002 FORMAT ( 14X24H VOLTAGE STEADY STATE = F13.8) 

1003 FORMAT (14X22h POWER STEADY STATE = F13.8//////) 

1004 FORMAT (//14X10HHTR ON F0RF9.4»4H HRS) 

3000 FORMAT ( 14F8.3) 

2000 CONTINUE 
END 
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